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FORiiARD 

Satellites have provided mati with a vantage point from Which to 
peer deeper into space as well as to retrospectively and spynoptically 
view his home planet. The SKYLAB earth orbital satellite has been the 
first manned U.S. flight laboratory to be devoted principally to remote 
observations and experimentation in the unique environment of space. 

The project was conceived in the early 1960's as an extension of the 
Apollo program which during that decade emphasized space exploration. 
During the period of the 60 's, application studies were consummated 
from which a proposed program for the conduct of solar observations 
and biomedical experiments evolved. In early 1970, the Apollo 
Applications Program was redesignated the SKYLAB Program and at about 
the same time the scope of the project was enlarged to include earth 
observation experiments. Earth observations were to be conducted with a 
variety of sensors including photographic cameras, an infrared spectro- 
meter, a multispectral scanner (visible and infrared), and a microwave 
radiometer/ scatterometer and altimeter. These sensors together with 
onboard support equipment were designated the Earth Resources Experiment 
Package (EREP). A total of 146 proposed investigations were selected to 
comprise the EREP data-user program. This paper reports oh the conduct 
and results of one of those investigations, specifically a study in the 
utilization of EREP data for monitoring changes in the breeding habitat 
of migratory Waterfowl. 

The authors wish to acknowledge Harvey K. Nelson who orginally 
conceived and developed the ideas behind this investigation. We are 
grateful too for the assistance of several individuals in the conduct 
of this work. In particular, we acknowledge the cooperation of our 
NASA technical monitor, Rigdon Joosten, the administrative guidance 
of VI. Reid Goforth and the field assistance of A. T. Klett, the 
latter two individuals both of the Northern Prairie Wildlife Research 
Center. Finally, but far from last, we are indebted for the efforts of 
Diana Rebel of the Environmental Research Institute of Michigan for her 
labors with data processing and analysis. 
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CHAPTER 1 

. : suMi^RY AND ;coweLysiONs;^^^ :■ 

To explore proceduv'es which could enhance the capabilities of the 
U.S. Fish and Wildlife Service for monitoring the breeding habitat of • 
niigratovy waterfowl j an BValuation of the SKYLAB Earth Resou^^ : " ^ 
Experiment Package (EREP) has been conducted. A related study had 
previously been carried out utilizing data collected by the LANDSAT-l 
satellite. The fact that the two studies have overlapped both ; ; 
chronologically and geographi cal ly has allowed the results of one to ■ 
reinforce the other and has al lowed for a comparison of the two sensor 
systems. In particular, m have emphasized the use of data collected 
by multi spectral scanners and the processing of these data using 
general purpose and special purpose digital computers. The use of 
automatic data processing techniques is uniquely suited to this type of 
task because of the wide expanse of prime waterfowl breeding areas and 
because of the need to quickly assimul ate and collate parametric 
information on habitat conditions. 

The specific objectives of both the LANDSAT and this SKYLAB study 
were to map and tabulate statistics on surface water conditions and 
to determi ne changes in wetness between the spring breeding period and 
the fledging period of July or early August in a glaciated prairie 
region in east-central North Dakota. This study has principally 
considered habitat conditions related to open surface water (i.e., 
ponds and lakes). The study as originally planned envisioned the use 
of two sets of data collected by the SKYLAB multi spectral scanner in 
May and again in July or early August. Because of operational constraints 
the SKYLAB/EREP System was not able to achieve this repetitive coverage. 
Instead a single observation occurred on 12 June 1973. The timing of 
the observation has allov/ed us to supplement the single SKYLAB obser- 
vation with bracketing LANDSAT observations which did occur on 14 May 
and 7 July of that year. This series of three sequential observations 
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provided an opportunity for tracking the dynamic water conditTons 
over nearly a two month interval. Approximately the same areas were 
observed on the three different dates. The study area included portions 
of two different physiographic regions, a Goteau or moraine feature 
created by; stagnation ice, and a drift plain or Tow relief feature of 
numerous ground moraines^ Because of differences in wetland frequency 
between them, these physiographic features have served as a basis for 
stratifying the numerical results obtained in this study. 

The mapping of open water has been carried out by two radically . 
different techniques , a si ngle-channel approach and a muTti pTe-channel 
approach termed "proportion estimation". The single channel approach 
delineates water by thresholdihg or level slicing the scale of radiance 
values in a single near-infrared waveband. This approach is effective 
because the radiance of water relative to other terrain features was 
uniformly low. Proportion estimation processing involved the use of 
multiple data channels and a computational algorithm for estimating 
the fractions of pure materials present within the resolution cell of 
a multi spectral scanner. , 

For the water delineation a single channel of near- infrared data 
was used to produce a computer generated thematic map and ralated 
statistics for an area encompassing 3833 km (1480 square miles). 

This area included portions of both the Goteau and Drift Plain physio- 
graphic regions. The results were produced with a single near-infrared 
waveband of data gathered by the SKY LAB multi spectral scanner (instrument 
experiment S-192). This multi spectral scanner included a total of five 
near-infrared wavebands in the 0,78-to 1 .75-um range, all of which had 
good signal-to-noise characteristics. Our experience has indicated that 
any one of these bands would have been reasonably satisfactory for small 
scale v/ater mapping by satellite. Given a choice, however, the 1.55- 
to 1.75-ym band was less ambiguous for water recognition and was easier 
to utilize in terms of training the computer. Use of the 1.55- to 
1.73-ym band also held the promise that shallow water features when re- 
solvable in the small scale satellite data Were more likely to be recog- 
nized as open surface water. 
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Although the maps which were generated in this effort graphically 
portrayed water conditi ons , without further interpretation and quant- 
ification, such displays do not fully satisfy management and research 
needs. Quantification- of the data was easily achieved when we used a 
digital computer to tabulate the size and location of each recognized 
water feature and thence to summarize the results. The summary indicated 
that the areal density of ponds in nearly a1 1 size classes considered 
(1 thru 50 acres) was generally an order of magnitude greater in the 
Coteau stratum than in the Drift Plain stratum. The data also indicate J 
that the SKYLAB enumeration of ponds as a whole was consistent with 
enumerations made using lANDSAT data. This is to say, that a decline 
in pond numbers noted from May to July indicated by the LANDSAT ob- 
servations was also reflected in the intermediate SKYLAB data. When 
individual ponds. as represented in the SKYLAB data were studied, however, 
we saw that their areal display did not follov; as consistent a pattern. 
Approximately 30 percent of the SKYLAB lakes examined were 1 arger than 
expected while another 30 percent were smaller than expected. We have 
concluded that these deviations from the mean tended to balance each 
other and that the source of the problem was the conical scanning con- 
figuration of the ..ensor and the procedures used to subsequently convert 
the data into a rectangular grid of scene elements or pixels. Both the,, 
scanning format and the associated techniques for data conversion 
appear to have had the net effect of slightly but systematically 
altering areal measurements and the geometric fidelity of small scale 
scene features such as prairie ponds. 

In another phase of the study, the limited testing of a unique 
technique for improving the apparent spatial resolution of multi spectral 
data was undertaken. The technique, termed “proportion estimation," 
involved the use of a computational algorithm for. estimating the fractions 
of pure materials present within the resol ution cell of a mul ti spectral 
scanner. To be effective, proporti on estimati on processing requires a 
high degree of spatial registration between spectral data channel s . Our 
early work with S KYLAB mul ti spectra 1 . scanner data in a li ne-straightened 
format (after conversion from a conical scan format) indicated the 
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occurrence of frequent jaut random channel to channel r..:sregistrat1on 
of as much as one or several pixels. The bulk of these registration 
errors were apparently introduced in the process of converting data 
froin a conical -scan to a line-straightened format. Consequently all 
of our mul ti spectral processing was accomplished with data in a coni cal - 
scan format although these data proved more difficult to handle and 
display. 

Proportion estimation computations were applied to an area of 286 

o 

km (no square miles). The resultant computer output was a set of 
water proportions for each scene pixel. Results obtained from pro- 
portion estimation processing indicated that the minimum discernable 
pond size was four-te minimum size detected with the single- 

waveband thresholding algorithm. Ponds smaller than this could be 
detected, however, numerous errors of false recognition (commission 
errors) also occurred when the finer degree spatial resolution was 
attempted. 

processing of LANDSAT data had previously 
2 

been applied to the same 286 km site^ It should be noted that nearly 
the same number of water features were identified in the LANDSAT data 
of 7 July 1973 as were observed with the SKYLAB observation data 
collected 25 days earlier on 12 June. If anything, this would indicate 
that the SKYLAB data did not yield as great a count of ponds as should 
be expected because some ponds should have diminished in size to the 
point of extinction over the interval. However, because of the 
deficiencies inherent in data which have undergone a format conversion, 
any inconsistencies of the small magnitude experienced could reasonably 
be attributed to the problems associated with data handling and 
conversion. 

It is most significant to note that certain lakes, which had been 
only partially delimited in the LANDSAT proportion estimation processing, 
were in fact fully delimited with the SKYLAB data. The lakes in 
question were shallow alkaline lakes ^w were high in suspended 
solids and/or precipitated alkali bottom sediments and which had an 
anomal ous appearance when compared to fresher lakes. We attri buted 
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this improved recognition capability to the availability and utilization 
of several wavebands in the near-infrared spectral region. The SKYLAB 
scanner provided six wavebands of near infrared information to a 
maximum wavelength of 2.35 urn whereas only two near-infrared wavebands 
to a maximum wavelength of 1.1 pm were available from the LANDSAT 
system. in particulars we feel that a waveband in the 1.5-to l.S-ym 
atmospheric window is important for the delineation of water and hygric 
scene features in general. 

Finally, we should point out that it is difficult to state con- 
clusively the results of this investigation because of numerous problems 
which became manifest during its course. Specifically, the two EREP 
observations upon which the measurement of habitat change was predicated 
were never realized. The single observation, which was obtained, did 
not occur during either the requested May or the requested July time 
periods. As a result the observation did not occur during an appropriate 
phonological period, and it did not coincide with supporting aircraft 
observations nor with the respective May and July breeding and production 
surveys routinely conducted by the U.S. Fish and Wildlife Service (FWS). 
These deviations from the originally planned experiment make it 
unrealistic and impractical for us to attempt at this time to assess 
the cost effectiveness of EREP data, as opposed to current techniques 
for predicting waterfowl populations. 

The diurnal timing of the observation was less than optimal in 
that it took place at 06:19 local solar time, too early in the day for 
sufficient illumination. As a result, the poor signal-to-noise ratios 
in the visible and thermal infrared wavebands precluded the plenary use 
of those spectral channels. Furthermore, the multispectral scanner 
observation only marginally encompassed the study site and many areas of 
environmental interest were not surveyed nor did the SKYLAB observation 
spatially coincide with many of the supporting aircraft and ground 
survey transects. These were problems of an operational nature 
occasioned by satellite system complications which arose during and 
after launch and by the need to accommodate the requirements of both 
this and many other investigations in a relatively short time. One 
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source of difficulty which was inherent to the scanner system and 
which has already been referred to was associated with the conical 
scanning confi guration and the subseguent data ref ormatti ng . 
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CONCLUSIONS 

1., Satemts remote sensing techniques hold Gonsiderable promise 
for the rapid, synoptic assessment of waterfowl breeding 
habitat. 

2. A simply implemented technique, requiring a single near- 
infrared waveband of data, exists, for delimiting open surface 
water features. This capability with the use of automatic 
data processing techniques has the potential for being 
operationally incorporated into ongoing habitat assessment 
programs in the near future. 

3. Trends noted in numbers and size distributions of water 
features were consistent between LANDSAT and. SKYLAB data sets 
when considering a large observation scene as a whole. 

4. The sizes of small ■individual water features were not, 
however, consistently determined. These areal inconsistencies 
of being rendered either too large or too small appeared to be 
due to the algorithm used to convert data from a curvilinear 
scan-line format to a straightened scan-line format and tended 
to average out when a 1 arge group of water features as a whole 
were considered. 

5. The proportion estimation technique, utilizing the added 
information content of multiple spectral wavebands, has alTowed 
for the recognition of a greater number of small ponds not 
previousl y i denti f i ed and al so great! y improved the area and 
peripheral shape definition of the larger ponds and lakes. 

The use of this technique has improved resolution capabilities 
for mapping open surface water by a factor Of between two to 
three over the nominal resolution limit of the data. Further 
testing and refinement of the technique will be requiredr 

6. Spectral information in the reflective infra'red (0,7 to 3.0 ym) 
is often unique and useful for the thematic classification of 
terrain features. In particular, we feel that a waveband in 
the 1.5-to 1.8-um «:mospher1c window is important for the 
delineation of water and hygric features in general. 
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The conical scanner while offering certain optical “mechanical 
advantages has in this program. created problems of a data 
handling and display nature. It has been necessary to convert 
the data from a curvilinear scan format into a rectangular or 
straightened scan-Tine format. This procedure has diminished 
the spatial accuracy and geometric fidelity of the data. 


2 

IMVESTIGATION BACKGROUND 

The primary breeding areas of North American waterfowl (primarily 
ducks) are the Dakotas, the southern portions of the prairie provinces, 
northwestern Canada and parts of Alaska (Figure 1). These areas of the 
mid-continent are major contributors in sustaining the total continental 
duck population which amounts upwards to 120 million birds. Small ponds 
and lakes of the glaciated prairie region, commonly referred to as 
prairie potholes, are the backbone of duck production habitat in North 
America (Figure 2). The prairie pothole region, composed of the 
southern portions of Alberta, Saskatchewan, and Manitoba and parts of 
North and South Dakota, Minnesota, and Montana makes up only 10 percent 
of the total breeding area of North America yet produces 50 percent of 
the continental duck crop in an average year (Smith, et al., 1964). 
Because of the area's Importance for waterfowl production the region is 
monitored annually by intensive systematic surveys conducted by the U.S. 
Fish and Wildlife Service in cooperation with the Canadian Wildlife 
Service and various states and provinces. Aerial surveys made in May 
and July are used with air-ground correction factors, to provide indices 
of breeding population size, habitat conditions, and waterfowl pro- 
duction. These indices serve to aid in making management decisions 
relating to annual hunting regulations and to answer certain research 
needs. Breeding ground survey data must be collected and summarized 
before early August when various national and regional waterfowl 
meetings convene to formulate annual hunting regulations. In some 
years breeding ground survey biologists are hard pressed to complete 
this task by early August. Descriptions of the operational aspects 
of the breeding ground surveys are given by Crissey (1957), Stewart, 
et al. (1958), and, more recently, Henny, et al. (1972). The use of 
survey statistics for modeling waterfowl production is discussed by 
Cooch (1969), Crissey (1969) and Geis, et al. (1969). 
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FIGURE 2. THE PRAIRIE POTHOLE REGION OF MID-CONTINENT NORTH AMERICA. 
(After Stewart and Kantrud, 1973.) The state of North Dakota is 

indicated by shading. 









Canadian and U.S. investigators have suggested the importance of 
pond numbers in regulating- annual waterfowl production (Cooch 1969, 
Crissey 1969). From the work of these and other researchers it has 
become evident that the degree of wetness, especially changes in pond 
numbers during May-July< should provide a useful index for predicting 
annual waterfowl production. In addition to annual estimates of pond 
numbers and distribution, information on short and long term ecological 
changes due to natural and economic causes are needed to properly 
manage this wildlife resource. A series of research investigations 
have and are being conducted on the premise that remote sensing 
techniques can be used to assimtlate such information on habitat 
conditions. 

A comprehensive program to investigate the_ potential applications 
of remote sensing techniques as a tool , in. the management and scientific 
study of waterfowl populations was developed in 1968 by the U.S. Fish 
and Wildlife Service and the Environmental Research Institute of 
Michigan (prior to 1972 known as the WITlow Run Laboratories of the 
University of Michigan). From 1968 to 1970 this work involved a 
series of airborne multispectral data collection and analysis ex- 
periments designed to assess waterfowl breeding habitat. The work 
was sponsored by the U.S. Department of the Interior's Earth Resources 
Observation Systems (EROS) Program and by the National Aeronautics 
and Space Administration (NASA). Biological interpretation and site 
coordination was provided by personnel of the U.S. Fish and Wildlife 
Service's Northern Prairie Wildlife Research Cen-ber at Jamestown, 

North Dakota. 

These experiments emphasized data collection by airborne 
multispectral scanners and data reduction and analysis with computer 
aided techniques specifically developed for this task. This early 
work was reported by Burge and Brown (1970), Nelson^ et al. (1970), 
and Work and Thomson (1974). A more recent program (Work 1974, and 
Work et al., in press) was conducted with LANDSAT-1 data. 

The LANDSAT study was a natural extension of the earlier aircraft 
program. Its goals were to: 1) modify and utilize techniques which 


evolved frotn the aircraft procjram, and 2) develop new techniques which 
would be suitable for high altitude, wide area Csynoptic) surveys. 

Data were to be used principally to document the amount of surface 
water present during a May and a successive duly observation. Data 
acquisition for the LANDSAT study was planned for the spring-summer 
of 1972. However, because of delays in the launch of the satellite, 
the bulk of that investigation was deferred until a May-July sequence 
of data was available from the 1973 season. As a result, the LANDSAT-1 
program and the SKYLAB investigation with which this report is concerned 
have utilized data observations which have been nearly coineidential 
in a temporal and spatial sense. Because of this unique situation we 
have made occasional references and comparisons in this report of 
data from both sensor systems. 



. . . DESCRIPTION OF THE INVESTIGATION . . . 

Our bveran goal for this inyestigatl on was to develop and test 
techniques for identifying and monitoring prairie waterfowl habitat 
using the EREP sensors. This and subsequent work is intended to bring 
us to our long term goal which is to define habitat quality indices 
which accurately relate to the dynamics of waterfowl production and 
which are discernabTe with remote sensing techniques. The application 
of high altitude (e.g. , spaceciraft) sensors and automated data handling 
is Uniquely suited to monitoring waterfowl habi tat because : (T) the 

survey is generally regional and international in scope, (2) the 
seasonal and annual changes in habitat conditions require repeated 
observations, and (3) the results are needed promptly for management 
decisions. 

In undertaking this study our first objective was to monitor 
changes in the breeding habitat of migratory Waterfowl between May, 
the peak nesting season for several species of ducks, and duly or 
early August when most duckl ing broods have hatched^ Proposed indi- 
cators of habitat quality were surface water, the general degree of 
terrain wetness, plant phenology, and land-use patterns. Primary 
emphasis, however, was placed upon the observations of ponds and 
lakes to include the assimilation of statistical data on their numbers, 
areal extent, frequency, distribution, and aggregation. Information 
on ponds and lakes is immediately important because it is currently 
used in models for predicting annual waterfowl production (^®isi et 
al. 1969). 

Consistent With earlier work (Burge and Brown 1970, and Work and 
Thomson 1974), the use of data obtained by a multi spectral scanner 
was stressed. Multispectral scanners offer the advantages of a multi- 
plicity of spectral wavebands and quantified data values in the form 
of digital signals recorded on magnetic; tape* The former advantage 
broadens the data's information content while the latter allows for 
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direct and rapid porcessing and analysis by digital computers. In 
addition to the primary role of scanner data in this investigation, 
photographi c imagery was uti 1 1 zed in a secondary capacity for the 
selection of computer training data and for the verification of recog*^ 
nition maps resulting from automatic data processing. 

The intended approach to carryi ng out this investi gati on was to 
use EREP multi spectral scanner data collected over the same site in 
a sequence of May and July observations. The resultant data were to 
be processed for the recognition of surface water features and for a 
determination of habitat change as evidenced by changes in surface 
water conditions. Expected output of the automatic data processing 
was thematic maps of ponds and a tabulati on of statistics on surface . 
water conditions. From the onset of the program it was recognized 
that the EREP scanner would not be capable of consistently delineating 
ponds less than about two hectares (5.0 acres). Therefore, in addition 
to merely mapping resolvable surface water features, estimates of the 
presence of ponds smal 1 er than the resol uti on limit were proposed by 
the use of a double sampling scheme in which low altitude aircraft 
data were to be collected at approximately the same time as each EREP 
data, pass. 

Because of operational delays in the launch and manning of the 
SKYLAB space station, data collection during the May 1973 breeding 
period did not occur. The subsequent compression of a large number of 
other pi anned earth observations i nto a foreshortened data acquisition 
period also precluded a second data observation as planned for July or 
early August of 1973. Instead a single data recording overflight 
occurred on 12 June 1973. Previously an overflight by supporting 
ai rcraft had occurred on 12 May 1973 and a subsequent aircraft overfl i ght 
took place on 12 August 1973 on schedules which coincided with both 
the waterfowl nesting and brood rearing seasons. We had hoped for 
near simultaneous SKYLAB. overflight and data observation. This 
failure to obtai n a repetition of seasonal observati ons and the un- 
timeliness of the single observation in terms of both season and 
coordination, with supporting aircraft necessitated a modification of 
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project goals and emphasis. Originally the program had been intended ; 
as a stand alone Investigation. Because of the circumstancesj however, 
it seemed appropriate to integrate this investigation with another 
ongoing study which was u ti 1 i zi ng LANDSAT-1 data. LANDSAT data 
: collected on 14 May 1973 and 7 July 1973 had previously been obtained 
and analyzed. It therefore was advantageous to utilize EREP data as . 
a third reference point in tracking the dynamic surface water conditions 
from May into July. 

■ Our 1 nability to obtain a temporal ly synchronized EREP and aircraft 
data set also precluded our conducting a double stage sampling experi- 
ment. In lieu of that experiment, we chose to utilize the manifold 
information of a multi spectral data-set to . systematically estimate 
percentage of surface v/ater present in each of a scene's resolution 
elements (pixels). This approach potentially resulted in the detection 
and tabulation of surface water features smaller than the sensor's 
nominal optical resolution limits. 

the text which follows includes a description of the. EREP multi- 
spectral scanner and a description of the study area and site data 
coverage (Chapter 4), foTlov/ed by the bulk of the technical report 
which discusses the methods and results of several analytic procedures 
used in this investigation. Chapter 5 is devoted to the mapping of 
surface water bodies; with a single vaveband of near- infrared data. 

This processing resulted in the generation of thematic water maps and 
statistics on water conditions and changes in these conditions. 

Chapter 6 is devoted to multi spectral data analysis particularly the 
use of multiple wavebands of data for estimating within-pixel surface 
;■ water content. Finally Chapter 1 contains a summary of investigation 
results and the conclusions derived from these results. 


- THE SKYLAB/EREP SURVEY OF WETLANDS 

m DAKOTA'. : 

This chapter provides additional background information regarding 
the principal EREP sensor system utilized for this investigation and 
characteristics of the study site. 



The SKYLAB/EREP Multispectral Scanner { Instrument Experiment S-1 92) 

■ was an optical '■mechani cal scanner coupl ed with a spectral dispersing 
and detector system/ The scanner assembly util ized a rotating .mirror 
scanning in the image plane of the collector optics to perform a conical 
scan of the object plane (i.e., the earth's surface). The cone angle 
was 5° 32/ about the instrument axis (nominally spacecraft nadir). The T 
spectral ly dispersed electromagnetic energy received from the earth ' s 
surface simultaneously irradiated 13 detectors, each detector being 
responsive to a unique spectral region. The scan pattern consisting 
of the forward 116°15' of the 360° scanning cycle covered a curvilinear 
path on the earth's surface with a swath or chord length of approximately 
72.4 km ( 39 .1 nautical mi les) and any desi red length along the ground 
track of the satellite (Figure 3), Approximately 94.8 scans occurred 
each second resulting in a scan-line to scan-line forward displacement : 
of approximately 72.4 m (238 ft,). . The sensor instantaneous-field-of- 
; view was 79.3 m (260 ft) square. Thus the instrument had an overscan 
of about 1 0%. Each , of the 13 detectors produced an electronic output 
signal corresponding to the average value of the radiance being 
received in its particular spectral band from a spot on the earth's 
surface contained in the instantaneous-field-of-view. The spectral 
range or waveband of each detector is given in Table 1. The analog ^ 
video signals enamating from each of the 13. spectral detectors were 
sampled and digitized at either high or low rates with the exception 
of the thermal infrared band (13) which was sampled at both rates. 

The low sampling rate corresponded to an approximate 72.6 m (238 ft.) 
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FIGURE 3. SKYLAB/EREU l-FULT I SPECTRAL SCANNER CONICAL SCAN PATTERN 


OP THE 

■ ‘-P 



TABLE 1 


SK.YLAB/EREP MULTISPECTRAL SCANNER (S-192) 
SPECTRAL RESPONSIVITY 


Band No. 

Description 

Nomina] Spectral 
Rangt (ym)* 

Measured Spectral 
Range (ym)** 

SDO 

Channel No. 

1 

Violet 

0.41-0.46 

0.420-0.447 

22 

2 

Violet-Blue 

0.46-0.51 

0.451-0.503 

18 

3 

Blue-Green 

0.52-0.56 

0.50-0.55 

1 & 2 

4 

Green-Yellow 

0.56-0.61 

0.54-0.60 

3 & 4 

5 

Orange-Red 

0.62-0.67 

0.599-0.654 

5 & 6 

6 

Red 

0.68-0.76 

0.654-0.734 

7 & 8 

7 

Infrared 

0.78-0.88 

0.770-0.890 

9 & 10 

8 

Infrared 

0.98-1.08 

0.930-1.050 

19 

9 

Infrared 

1.09-1.19 

1.030-1.190 

20 

10 

Infrared 

1.20-1.30 

1.150-1.280 

17 

11 

Infrared 

1.55-1.75 

1.550-1.730 

11 & 12 

12 

Infrared 

2.10-2.35 

2.10-2.34 

13 & 14 

13 

Thermal Infrared 10.2-12.5 


15,16 & 21 


" The nominal spectral range is referenced throughout this text. 

** S-192 spectral response calibration per National Aeronautics and 
Space Administration (1974). 
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center-to-center spacing and was used for bands 1, 2, 8, 9, 10, and 
13. The high sampling rate, equivalent to twice the low rate, approx- 
imated a 36.3 m (119 ft.) center-to-center sample spacing. Bands 3, 

4, '5, 6, 7, n, 12, and 13 were sampled at the high rate with even and 
odd numbered samples being handled as two low rate channels thereafter. 
The channels are referred to as SCO's (Scientific Data Output) in this 
report. In theory, all even numbered SCO’s were in spatial registration 
with each other. Similarly all odd numbered SCO's were in spatial 
registration but one-half pixel out of registration with the even SDO's. 

For the convenience of those investigators who utilized S-192 data 
and who desired data in a straight scan-line format, the conical data 
were available in a transformed line-straightened format. The algorithm 
used for this transformation was based on a nearest neighbor decision 
rule. After scan-line straightening, each high rate band was separated 
again into two low rate channels (SDO's) with the odd elements going 
into one SCO and the even elements going into another SCO as with the 
conical data described above. This investigation has utilized data in 
both a line-straightened format (Chapter 5) and a conical format 
(Chapter 6). 

SKYLAB/EREP Mission Profile 

The ^■.ission consisted of series of four launchings -- the 

first, to ir-^f ■ r the large unmanned laboratory in earth orbit and there- 
after three serial launches to carry crews to the orbiting laboratory. 
Each crew, consisting of three men transited from earth to the space 
laboratory, occupied the laboratory for an extended period, and sub- 
sequently returned to earth. The space laboratory was occupied for 
periods of 28, 60, and 85 days with intervals of 36 and 51 days respec- 
tively between the manned periods. This investigation utilized data 
collected during the first manned period which lasted from 25 May 1973 
to 22 June 1973. During each occupation, the crews conducted a series 
of biomedical, astronomical, engineering, and earth survey experiments 
all of which had to be closely managed to fit within a tight mission 
schedule. In addition to the necessary imposition of a time budget. 


other factors also posed restrictions to the earth resources survey 
program. They included spacecraft power limitations, the availability 
of consumables such as film and magnetic tape, spacecraft orbit location 
with respect to targets of interest, and the occurrence of suitable 
solar illumination and atmospheric viewing conditions. 

Throughout its useful life the SKYLAB Orbital V/orkshop circumvented 
the earth in a circular orbit inclined at approximately 50° to the 
equator. The orbit was controlled to yield a five-day repeating ground 
track. However, it should be noted that the five-day orbital re-traces 
progressed tnrough periods of darkness and daylight and that orbital 
drift, which at times became considerable, did occur. 

Typical orbital paths projected onto the earth’s surface are shown 
in Figure 4. The single earth observation made on behalf of this 
investigation occurred on 12 June 1973 on ar\ ascending node of the 
orbital track crossing the state of North Dakota from the northwest 
corner thence east-south-east passing approximately over Fargo. Table 
2 lists the parameters of this overpass. 

The North Dakota Study Area 

The study site chosen for this investigation, although small 
relative to the regional scope of the prairie pothole country was 
extensive enough to provide an adequate test of survey methods. The 
specified intensive study area was centered on Woodworth Station* and 
extended eastward to longitude 100°00*. The test site was situated 
completely within the North Dakota prairie pothole biotic area but 
did encompass two distinctly different groups of glacial landforms or 
physiographic divisions — the Missouri Coteau and a Glacial Drift 
Plain. Figure 5 shows the biotic areas of North Dakota and the location 
of Woodworth Station. The Coteau overlying approximately two-thirds 


*Woodworth Station is a field research site operated by the 
Northern Prairie Wildlife Research Center ofpthe U.S. Fish and Wildlife 
Service. The Station, approximately 15.5 km*^ (6 in size, lies 5 
km (3 miles) east of the village of VJoodworth, North Dakota or 48 km 
(30 miles) northwest of Jamestown. The station coordinates are 47°08‘N 
and 99°14'W. 


22 



REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


TABLE 2 

SKYLAB MULTI SPECTRAL SCANNER (S-192) CHARACTERISTICS 


Scan Format; 

Scan Cone Angle; 

Active Scan: 

Scanner Optical 
InstantaneouS“Field“of-Vi£iW: 

Scan Rate; 

Altitude: 

Ground Radius of Scan: 
Sampling Rate (Along Scan); 

Number of Samples per Scan; 

Ground Speed: 

Satellite Ground Distance 
Forward Per Scan; 

Analog to Digital Conversion: 


Conical 

5 “ 32 ’ 

Forward 116°15’ of Scan 

0.182 rarad (79.71 m @ Altitude Shox<m) 

94.792 Scans /Second 

437,957m* 

42427 m* 

Loxif Rate Channels — 

72.6 m Center to Center 
High Rate Channels — 

36.3 m Center to Center 

Law Rate Channels - 1240 
High Rate Channels - 2480 

6866 ra/sec* 

72.43 m* 

8 Bit Words 


* These parameters apply specifically to the North Dakota overpass 
of 12 June 1973. 
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FIGURE 5. THE BIOTIC REGIONS OF NORTH DAKOTA. (After Stewart and Kantrud, 1973) 
The location of Woodworth Station with respect to the Missouri Coteau and a drift 

plain feature is Indicated. 
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of the study area is characterized by the prominence of high underlying 
bedrock which acted to buttress the advance of Pleistocene ice sheets 
causing extensive glacial stagnation. In addition, the ice had become 
overlain with large quantities of superglacial till which caused 
protracted and irregular melting of the underlying ice. The resultant 
topography is hummocky, drainage is non-integrated, and lakes and 
sloughs are abundant as is characteristic of an area of collasped ice 
topography. The Coteau is said to contain some of the best waterfowl 
breeding areas in the 48 contiguous states (Clayton, 1967). The drift 
plain was formed by glaciation that possessed a margin which retreated 
in an orderly manner and which occasionally halted or readvanced. 
Drainage in the plain is integrated only along the edges of the large 
melt water channels. Numerous shallow, marshy depressions are present 
between these former channels. Relative to the Coteau, however, the 
Drift Plain has fewer potholes and because of its low relief has been 
subjected to numerous wetland drainage projects. This difference in 
wetland occurrence warranted a stratification of the statistical results 
in this investigation. 

The areal extent of the multi spectral scanner data utilized in 
this study is shown in Figure 6. Generally data which were common to 
three observations (two by LANDSAT and one by SKYLAB/EREP) are presented 
in this report. For the EREP S-192 observation, the spacecraft's 
ground track passed within 34 km (21 statute miles) of Woodworth 
Station on a heading approximately 109°. Because of the relatively 
narrow lateral field-of-view of the EREP multi spectral scanner and the 
short data-take period, Woodworth Station was on the extreme margin of 
the scanner's coverage and scanner coverage only minimally included 
areas lying within the Coteau physiographic division which lay generally 
to the west and southwest of Woodworth Station. It is unfortunate that 
many of the sampling transects flown by the supporting aircraft were 
also outside the area actually observed by the EREP S-192 scanner. The 
short along-track duration of the S-192 observation, however, was 
necessitated by the limited supply of magnetic recording tape carried 
aboard the spacecraft. 
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FIGURE 6. GEOGRAPHIC LOCATION OF SATELLITE AND AIRCRAFT OBSERVATIONS. The larger 
lakes and major streamcourses within the state of North Dakota are indicated. 
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detection OF SURFACE WATER FEA 

SINGLE WAVEBAND OF NEAR-INFRARED' DATA 

LeveT thresholding of a radiation signal in a single near-infrared 
waveband is a reliable and simple technique for delineating surface 
water. This technique is effective because at hear-ihfrared wavelengths 
the apparent radiation of water is usually uniform and lower than for 
other terrain objects i. Thus using an appropriate hear-infrared waveband j 
water may be delineated by accepting scene points with low radiance 
values (classified as water) while rejecting all values above a certain 
threshold (non-water). 

In this texts we have termed this form of surface water detection 
"thresholding". To appreciate the effectiveness of thresholding, it 
is helpful to have an understanding of the interaction of incident 
radiation* with water. 

Spectral Radiance of Water 

The apparent radiance of a body of water is the result of: (1) 

reflections at the air-water interface, (2) reflections from particulate 
matter suspended in the water, and (3) reflections from the bottom. 
Because the fields-of-view of the LANDSAT and SKYLAB multi spectral 
scanners have been limited to near vertical observations and because 
water surfaces reflect specularly, radiation reflected by water to the 
scanners must have emanated from a sky position near the zenith. Given 
the northerly latitudes which characterize the glaciated prairies, 
satellite scanners viewing only near nadir generally do not view water- 
reflected direct solar radiation (i.e. the ground specular point is a 
considerable distance outside the field-of-view of the scanner). This 
leaves only that fraction of diffuse skylight which emanates from a 
near-zenith sky location to impinge upon the water surface and thence 


*This discussion excludes consideration of thermal or self-emitted 
infrared radiation and is therefore limited to radiation in the visible 
and near-infrared at wavelengths somewhat less than 4 ym. 
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to bo reflected to the scanner. In relative magnitude; however, 
diffuse skylight is much weaker than direct solar radiation, especially 
in the near-infrared and under clear sky conditions when optimal 
satellite observations are ppssible.: McDowell (1974) illustrated 
the magnitude and spectral differences between diffuse skylight apd , 
direct solar radiation (Figure 7). In addition, the reflected skylight 
component is further diminished because V'/ater surfaces are a uniformly 
weak ref lector of radiation which impinges at any but very oblique . 
angles (Figure 8), 

in considering reflections emanating from particulates within the 
water volume and from the bottom surface, water's absorptivity must be 
considered. In the near-infrared, that fraction of radiation which 
penetrates the air-water interface is largely absorbed, the extent of 
absorption being dependent upon the wavelength and the length of the 
water path. This situation is shown quantitatively in Figure 9 which 
illustrates the spectral transmission of pure water for a variety of 
path lengths. Consequently, a sensor viewing a water body in a near- 
infrared band receives little or no radiation that may have been re- 
flected by the bottom, or volume suspended particulates. 

In an earlier study utilizing aircraft data. Work and Thomson 
(1974) evaluated the relative merltS; of various near infrared bands 
for mapping surface water- They compared bands in the 0.73- to 0,92-ym, 
1.0- to 1.4-ym, and 1.5- to 1.8-ym ranges and concluded that all pro- 
duced reasonably good results. Given a choice, however, the longer 
wavelength bands did provide some marginal improvement. Longer wave- 
length alone should not predicate the choice of a water mapping wave- 
band however. For example, the use of a waveband in the 2.0- to 
2.6-]jm atmospheric window is not optimal because of the decreasing 
amount of solar radiation at. these wavelengths. It must be remembered 
that most terrestrial objects are relatively strong diffuse reflectors 
and that with adequate solar illumination such targets will contrast 
sharply with surface water features which consistently are darker. Work 
and Thomson (1974) concluded that an ideal waveband for delineating 
surface water lay within the 1.5- to 1.8-ym atmospheric window. 
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J 10 


Direct Solar Light + Skylight 
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FIGURE 7. SCENE lERADIANCE COMPONENTS FOR A CLEAR DAY. (After McDowell, 
1974.) Spectral irradiance levels attributable to skylight alone may be 
significantly different and variable with wavelength due to variations in 
atmospheric haze. The conditions shotm were recorded on an exceptionally 
clear day in New Mexico on 30 October 1970 at a sun elevation angle of 43 
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ANGLE OF INCIDENCE MEASURED 
FROM NORMAL 


FIGURE a. VARIATION OF REFLECTING POWER OF THE AIR-WATER 
INTERFACE AS A FUNCTION OF THE ANGLE OP INCIDENCE (Natural 
or unpolarized light of wavelength 0.589 pm) , 
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FIGDEE 9. SPEGTEAL TRAl^ISMITTANGE OF PURE WATER FOR DIFFERENT PATH 
LENGTHS. (Plotted after data from Sverdrup et al., 1942.) 


A graphic in ustratl on bf the- relative utility of different 
near-infrared bands for detecting water may be seen in Figure 10. 
Illustrated is imagery acquired by an airborne multispectral scanner 
flown in support of this study. For the three infrared bands illustrated 
(0.67- to 0.94-vm, and 1.0- to 1.4-ijm), standing water is rendered as a 
dark feature in the imagery due to the low level of radiant energy 
emanating therefrom. These wavebands are generally comparable to the 
SKYLAB multispectral scanner bands 7, 10, and 11 respectively (0.78- to 
0.89-ym, 1.20- to 1.30-pm, and 1.55 to 1.75-nm). 

The 0.67- to 0.94-ym imagery of Figure 10 illustrates an anomalous 
condition for at least one pond. The pond labeled "a" contained several 
light-toned, pincer-shaped features which occurred within the pond 
perimeter but which did not appear in the imagery at wavelengths greater 
than 1.0 mm. These features were due to a floating algal mat and 
possibly to plant submergents (water-milfoil, Myriophyllum exalbescens , 
and bladderwort, Utricularia vulgaris ! which may have been exposed by 
low water levels. In the 0.67- to 0.94-pm band, this vegetation was a 
moderately strong reflector of incident radiation, and consequently a 
light-toned rendition occurred in the imagery thus masking the under- 
lying water. 

In general, the radiance of vegetation is largely affected by the 
critical reflection of incident light from cell walls within a leaf 
(Gates, et al,, 1965, Gausman, 1974). Although radiation may be 
reflected several times before leaving the leaf, most of the radiation 
will be returned if there is little absorption by leaf tissues. This 
would seem to be true for the floating vegetation in the 0.67- to 
0.96-iim waveband of Figure 10. However, beyond 1.0 nm, water's increasing 
near infrared absorptance (per Figure 9} appears to influence the 
radiance of plant materials. In the 1.0- to 1.4-ym imagery of Figure 
10, the floating vegetation was no longer discernable from its water 
background due, possibly, to the high moisture content of aquatic 
plant tissues. In addition, the fact that the algal mat may have been 
floating several millimeters below the pond surface would also have 
precluded a string return of radiation. 
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Bare Soil "a" -Floating Vegetation 
Video Reproduction of 0.67-0.94fim Scanner Data 


Bare Soil 

Video Reproduction of 1.0-1.4/xm Scanner Data 



"c"-Shelterbelt "b" -Marsh 

Video Reproduction of 1.5-1.8/im Scanner Data 


-Shelterbelt 



FIGURE 10. COMPARISON OF SEVERAL NEAR-INFRARED WAVEBANDS FOR RENDERING 
STANDING WATER AND OTHER MOISTURE RELATED CONDITIONS. The floating vegetation 
at location ”a" partly obscured the standing water in the 0.67-0.94/im data only. The marshes 
at locations "b" and the shelte. belt trees at locations "c" had foliar vegetation containing liquid 
water. Many of the upland herbaceous plants were highly desiccated. These hygric to .xeric 
conditions are discernible by tonal differences in the 1,5-1.8/im data. The above data were 
collected by an airborne multispectral scanner operated at an altitude of 4500ft. above the Wood- 

worth Station on 12 August 1973, 1633 GMT. 
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Liquid water within a leaf is largely the cause of moderate ab- 
sorption beyond 1.4 um and according to Gates, et al . (1965) very 
strong absorption beyond 2.0 urn. Other investigators including Olson 
(1969), Rohde and Olson (1970), Myers, et al . (1970), and Gausman 
(1974) have demonstrated the influence of leaf water in suppressing 
reflectivity beyond 1.5 }jm. This situation is illustrated at locations 
"b" and ”c" in the 1.5- to 1.8-tjm imagery of Figure 10. At locations 
labeled "b", marsh vegetation consisting predominantly of bulrushes 
( Scirpus spp .) and cattails ( Typha spp . ) was present. Deciduous trees 
comprised the shelterbelts at locations labeled "c". Both of these 
communities had received some moisture during the usually dry summer, 
either because they were deep rooted (i.e., the trees) or because they 
were rooted in standing water (i.e., the marsh vegetation). Gross 
differences in moisture conditions did exist in the scene. Specifically, 
an extremely desiccated condition is exhibited by much of the dry 
herbaceous upland vegetation, most of which was dead. Reflections 
from these materials were high, and they contrast sharply with the 
marsh and shelterbelt communities in the 1.5- to 1.8-um imagery. 

The water recognition map included in Figure 10 was generated by 
thresholding 1.5- to 1.8-um data. The recognition map illustrates that 
in spite of the apparent low radiance of water, marsh, and shelterbelt 
features, standing water was unique for its low radiometric signature. 

Had water been mapped by thresholding the 0.67- to 0.94-ijm band or any 
of the SKYLAB multispectral scanner bands between 0.7 and l.Oum, the 
pond at "a" and other similarly occluded or shallow water features would 
at best have been only partially recognized. 

Methods 

Implementation of the thresholding technique was accomplished by 
observing radiance ilues for known water features within a scene and 
comparing these values with those of other terrain features also known 
to exhibit relatively low radiance characteristics. A decision boundary 
or threshold was then- selected which effectively separated surface water 
from all other scene objects on the basis of their relative brightnesses 
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(i.e., differences in apparent radiance). Experience has shown that 
various scene objects exhibit radiances that may or may not be close to 
the low radiance of water depending on (1) the specific near-infrared 
band under consideration, (2) the geographic locale, and (3) the 
phenol ogic state of some scene objects. 

In eastern North Dakota, dark prairie soils, Mollisols {formerly 
referred to as Chernozems), have consistently approached the low 
radiance values of water in near-infrared wavebands of less than 1.4 urn. 
In the 1.5- to 1.8-pm atmospheric window, vigorous green vegetation 
most nearly approach water's low apparent radiance (these water, soil, 
and green vegetation radiance characteristics were previously illustrated 
in Figure 10). For our North Dakota study site, the practice was to 
threshold or differentiate between water and either bare soil or 
vigorous green vegetation depending upon the particular near-infrared 
waveband utilized. 

The SKYLAB multi spectral scanner had five wavebands in the wave- 
length range 0.78 to 1.75 ym, any of which were potentially useful for 
discriminating open surface water using the thresholding technique. 

Thus the SKYLAB data offered a further opportunity to appraise the rela- 
tive usefulness of several near-infrared wavebands. Such an evaluation 
was conducted, the results of which are presented in Figures 11 through 
14. In the waveband considered in Figures 11 through 13 (0.78- to 
0,88-pm, 0.98- to 1.08-um, and 1.20- to 1.30-um respectively) , bare 
soil was the terrain feature most likely to be mistaken for open surface 
water. Each histogram represents a sample size of about 350 pixels 
drawn from throughout the observation scene. The overlap of water and 
bare soil histogram tails generally decreases with increasing wavelength. 
In the 1.55- to 1.75-pm waveband (Figure 14), the terrain material most 
likely to be mistaken for open water was vigorous green vegetation; but 
there was no overlap of histograms in this waveband. Our conclusion, re- 
affirmed by these data, was that the 1.55- to 1.75-ym channel was the 
least ambiguous for water discrimination. For this waveband a threshold 
boundary of 10 volts and less was selected for the delineation of open 
surface water. In affixing units to this threshold value, no inference 
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20 - 



FIGURE 11. HISTOGRAMS OF REFLECTANCE VALUES IN THE 0.78- to 0.88- pm WAVE- 
BAND. The distributions each represent approximately 350 pixels sampled 
from throughout a SKYLAB multispectral scanner observation made 12 June 

1973 over eastern North Dakota. 


20 - 

18- 



FIGURE 12. HISTOGRAMS OF REFLECTANCE VALUES IN THE 0,98- to 1.08- pm WAVE- 
BAND. The distributions each represent approximately 350 pixels sampled 
from throughout a SKYLAB multispectral scanner observation made 12 June 

1973 over eastern North Dakota. 
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FIGURE 13. HISTOGRAMS OF REFLECTANCE VALUES IN THE 1.20- to 1.30- pm WAVE- 
BAND. The distributions each represent approximately 350 pixels sampled 
from throughout a Slcylab multispectral scanner observation made 12 June 

1973 over eastern Worth Dakota. 
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10 20 30 

DIGITAL (VOLTAGE) COUNT - SDO 11 (1.55-1.75 pm) 

FIGURE 14. HISTOGRAMS OF REFLECTANCE VALUES IN THE 1.55- to 1.75-pm WAVE- 
BAND. The distributions each represent approximately 350 pixels sampled 
from throughout a Skylab multispectral scanner observation made 12 June 

1973 over eastern North Dakota. 
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of radiance values was intended. Indeed, we have only been concerned 
with relative differences of radiance between target types and con- 
sequently have made no attempt, using available calibration data, to 
relate the digital (voltage) count recorded on the data tapes to 
absolute radiance levels. i 

In this investigation, the selection of an appropriate threshold | 

boundary was manually accomplished, based upon a visual examination of 
a computer generated statistical analysis {i.e., histograms or frequency 
distributions). Given the limited amount of data to be analyzed and j 

the clear separability between terrain classes (i.e., open surface water \ 

versus vigorous green vegetation), this manual interaction in the | 

recognition process has been sati factory. The procedure could have j 

been automated, however, and the computer could have been used to | 

perform a one channel linear recognition between two object classes [ 

in a manner similar to maximum likelihood classification recognition | 

(or pattern recognition). Such an approach should be appropriate when 
it is necessary to handle large amounts of data which represent a 
variety of scene and illumination conditions. 

Results of Single Channel Water Recognition 

A computer-generated thematic map identifying open surface water 
over a 3621 km^ (1397 mi^) area as observed by the SKYLAB/EREP scanner 
was produced by thresholding a 1.55- to 1.75-Mm waveband (Figure 15). 

The tract shown partially overlapped both the Missouri Coteau (257o) 
and the Drift Plain (75%) strata. It is evident from this map that the 
Coteau (lower left) had a considerably higher density of ponds and 
lakes. The larger lakes in the Drift Plain were frequently the result 
of major imf 'undments on the James River watercourse (Jamestown 
Reservoir, Jim Lake, Mud Lake, Arrowwood Lake, and Juanita Lake). 

The SKYLAB/EREP observation of 12 June 1973 was interleaved between 
two LANDSAT observations which occurred on 14 May 1973 and 7 July 1973. 

Thematic water maps resulting from those LANDSAT observations are shown 
in Figures 16 and 17 respectively. Only the upper 40 percent of each of 
the LANDSAT maps is common to the area mapped by SKYLAB. Note, too, 
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that the maps scales differ. The LANDSAT maps were Initally generated 
at a scale of l/24s000 using a general purpose, digital printer writing 
132 pixel columns per each 15-inch wide paper strip. The LANDSAT maps 
illustrated each consisted of six printed strips which were manually 
abutted and subsequently photo-reduced to a convenient size. The SKYLAB 
map was initally generated at a scale of 1/155,000 using a computer 
controlled, ink-jet printer. This format was of a more convenient size 
and was considerably easier to reproduce. Thus the ink-jet map repre- 
sented a significant reduction in production time and labor. 

The three sequential observations represented by Figures 15, 16, 
and 17 spanned a time period of low precipitation and progressive 
desiccation. The diminishing of prairie ponds and lakes is normally 
to be expected during a May to July period. Visually, however, it is 
difficult to discern from the maps what the changes in surface water 
conditions were, and without extensive manual interpretation it is 
impossible to quantify these changes. For purposes of analysis, 
therefore, the data are more conveniently assessed if they are assembled 
in statistical form by automatic data processing techniques. Figure 
18 illustrates such a statistical tabulation as generated from tape 
recorded data gathered by the SKYLAB scanner over the Coteau physio- 
graphic stratum. The upper tabulation (only partially illustrated) 
is an enumeration of all recognized surface water features while the 
lower tabulation summarizes the frequency of pond occurrence by size. 

For the biologist concerned with the management of waterfowl populations 
this summary provides a ready assessment of habitat conditions over 
wide (synoptic) areas. The SKYLAB map and the related statistics have 
been generated from tape recorded data which was made available to us 
in a scan-line straightened format. Existing software and printers 
were capable of handling data in only this format. 

The statistical tabulations were produced by a software program 
designed specifically for use with satellite data. The program was 
adapted from a pre-existing program which was used with wide field-of- 
view, low altitude (aircraft) scanner data. These programs function 
by the use of a threshold decision criterion for classifying a grouping 
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FIGURE 18. EXAMPLE OF COMPUTER PRINTOUT OF POND AND LAKE STATISTICS FOR 
AN AREA WITHIN THE MISSOURI COTEAU PHYSIOGRAPHIC PROVINCE OF NORTH DAKOTA 
The above statistics resulted from a SKYLAB multlspectral scanner (S-192) 

observation of 12 June 1973. 




of pixels as a water body arid thereafter computing the area of each 
specific water feature. Perimeter and shape factor (a measure of 
shape complexity} cal cuTations were also possible. Because of the 
small size of many of the ponds and the potential that the shorelines 
of larger ponds and lakes could vary widely In length at scales finer 
than the resolution limit of the data, perimeter and shape factor 
calculations were not performed. Appendix A of this report describes 
in detail the computer program used to generate these statlcticM^^^^ " 
data while Appendix B includes a complete tabulation of surface 
water statistics derived from the SKYLA8 scanner observation of 12 June 
1973. 

In the statistical tabulations, we Identified each pond and lake 
and defined its position with each of two coordinate systems based on 
a scan-line and point number scheme and a more conventional latitude 
and longitude system. The scan-line and point number information was 
inherently available from the digital tapes. Conversion to latitude 
and longitude coordinates was accompTi shed by a regression analysis 
which used several control points located within the scene. The 
convention used to reference each water body was to identify the body 
by the number of the last scan-line with at least one pixel in the 
water body and the point number of the greatest numbered water pixel 
of that scan “line. 

Graphical summaries of pond frequency for the two observed strata 
are shown in Figure 19. These data were normalized so that comparisons 
between the different sized strata are possible. Note that the ordinate 
or density scales for the two strata are logarithmic and that the two 
density distributions differ between strata by approximately an order of 
magnitude in all size classes. A summary of seasonal change in pond 
numbers over several consecutive time intervals is illustrated in 
Figures 20 and 21. These changes were observed in the Coteau stratum 
by the use of a combination of LANDSAT and SKYLAB observations*. The 


*A similar measure of surface water change for the Drift Plain 
is not available because the processed LANDSAT data did not represent 
a large enough sample in this stratum. 


DENSITY (NUMBER OF PONDS PER SQUARE MILE) 



FIGURE 19. SUMMARY OF SIZE DISTRIBUTION OF PONDS IN THE GOTEAU AND 
DRIFT PLAIN STRATA AS DETERMINED USING SKYLAB MULTISPEGTRAL SGANNER 
DATA GOLLECTED 12 JUNE 1973. I-Ihere the pond size Increments are 
greater than one acre, the data have been averaged over the increment. 
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APPARENT POND SIZE (ACRES) 

FIGURE 21. SUM^^ARY OF SIZE DISTRIBUTION OF PONDS IN THE COTEAU STRATUM FOR A 25-DAY 
LATE SPRING/SUMMER INTERVAL. Where the pond size increments are greater than one 
acre, the data have been averaged over the increment. 
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graphical data indicate a progressive drying trend for May thru June 
to July. However, the trend appeared to be protracted during the 14 
May to 12 June interval after which ponds diminished in size and 
numbers at a more rapid rate. On a synoptic basis, at least, it appears 
that the two sensor systems (i.e. , SKYLAB and LANDSAT) were consistent 
in that both data sets indicated a decline in both area and numbers of 
surface water features. It should be noted, however, that the data 
are for sampled areas which were not specifically the same for each 
of the three observations. As a result, the observed trends for the 
larger ponds and lakes were not always consistent because of the low 
sampling frequency. 

The trends, as noted above, appeared to be manifested in both 
the SKYLAB and L^^NnSAT observations when each of the data sets was 
considered as a whole. One may then ask whether individual ponds 
and lakes as observed by both the SKYLAB and LANDSAT scanner systems 
also adhered to the group trend. To answer this question, 21 ponds 
and lakes ranging in size from 2 to over 190 hectares were randomly 
selected for comparison. Figure 22 is an enlargement of a section 
of the thematic water map shown previously in Figure 15. Indicated 
in Figure 22 are the names and locations of the 21 ponds and lakes 
studied. Table 3 is a listing of these water bodies along with their 
geographic coordinates as listed in the computer output stream. Note 
that the coordinates are consistent between data sets but that slight 
differences exist between the SKYLAB and LANDSAT observations. These 
discrepancies can be attributed to the fact that the computed coordi- 
nates from the SKYLAB data were for the northeastern corner of each 
water body while the LANDSAT data were for the southeastern corner. 

This difference is particularly apparent for the larger water bodies, 
for example, Barnes Lake. 

Table 4 is a listing of water body size for the same 21 ponds 
and lakes. The LANDSAT data indicated a decline in pond size from May 
to July in all but two instances, Lawrence Lake and Fish Lake, which 
virtually remained constant. On the other hand, the SKYLAB pond and 
lake data which were collected intermediate between the LANDSAT 










TABLE 3. 

COMPARISON OF COMPUTED LATITUDE & LONGITUDE COORDINATES 
PER LANDSAT & SKYLAB OBSERVATIONS 



Computed Latitude (Deg. N) 
Per Observation of: 

Computed Longitude (Deg. W) 
Per Observation of: 


LANDSAT 

5/14/73 

LANDSAT 

7/7/73 

SKYLAB 

6/12/73 

LANDSAT 

5/14/73 

LANDSAT 

7/7/73 

SKYLAB 

6/12/7: 

Barnes Lake 

47.224 

47.225 

47.239 

99.286 

99.286 

99.260 

Norden Lake 

47.255 

47.255 

47.260 

99.258 

99.258 

99.255 

Jerome Lake 

47.246 

47.246 

47.248 

99.261 

99.261 

99.259 

Janice Lake 

47.242 

47.242 

47.246 

99.250 

99.250 

99.249 

Colby Lake 

47.208 

47.208 

47.211 

99.262 

99.263 

99.260 

Northwestern 

Lake 

47.210 

47.210 

47.214 

99.217 

99.215 

99.206 

Alkali Lake 

47.204 

47.205 

47.208 

99.205 

99.205 

99.203 

Trautman Lk. 

47.177 

47.178 

47,179 

99.200 

99.200 

99.201 

Hall Lake 

47.167 

47.168 

47.172 

99.250 

99.250 

99.237 

Hotchkiss 

Pond 

47.151 

47.151 

47.153 

99.263 

99.264 

99.262 

Clark Lake 

47,134 

47.135 

47,136 

99.238 

99.236 

99.236 

Fish Lake 

47.130 

47.130 

47.130 

99.234 

99.231 

99,229 

Big Lake 

47,135 

47.136 

47.138 

99.229 

99.227 

99.226 

Goldwin Lake 

47.136 

47.137 

47.140 

99.196 

99.196 

99.195 

Limesand 

Pond 

47 . 139 

47.140 

47.140 

99.182 

99.182 

99.182 

Schelske Lk. 

47.108 

47.111 

47.111 

99 . 196 

99.195 

99.194 

Bust Lake 

47.200 

47.201 

47.206 

99.293 

99.294 

99.289 

West Lake 

47.203 

47.204 

47.209 

99.232 

99.231 

99.297 

Lawrence Lk, 

47.208 

47.208 

47.212 

99.168 

99.169 

‘>9,168 

Eugene Pond 

47.175 

47.176 

47.177 

99.174 

99.174 

99.174 

Woodworth 

Marsh 

47.134 

47.135 

47.135 

99.300 

99.300 

99.302 
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TABLE 4. 

COMPARISON OF COMPUTED WATER AREAS PER 
LANDSAT & SKYLAB OBSERVATIONS 




(a) (b) (c) (d) (e) 



Computed Surface Water 
Areas (Hectares) Per 
Observation of: 

LANDSAT Observed 
Mean Area 

SKYLAB 

Departure 

from 

LANDSAT 

Mean 

(No. of 
Pixels) 

LANDSAT 

5/14/73 

SKYLAB 

6/12/73 

LANDSAT 

7/7/73 

(Hectares) 

(No. of 
Pixels) 

Barnes Lake 

191.15 

197.17 

180.01 

185.58 

368 

+18.24 

Norden Lake 

33.11 

35.76 

30.67 

31.89 

63 

+6.09 

Jerome Lake 

16.78 

17.38 

13.78 

15.28 

30 

+3.02 

Janice Lake 

17.66 

17.38 

16.45 

17.06 

34 

+0.52 

Colby Lake 

22.07 

23.84 

20.45 

21.26 

42 

+4.06 

Northwestern Lk. 

42.38 

38.74 

37.34 

39.86 

79 

"1.76 

Alkali Lake 

29.14 

26.32 

24.45 

26.80 

53 

-0.74 

Trautman Lake 

8.83 

3.97 

6.22 

7.53 

15 

-5.59 

Hall Lake 

32.67 

27.32 

23.11 

27.89 

55 

-0.76 

Hotchkiss Pond 

7.06 

3.48 

4.00 

5.53 

11 

-3.23 

Clark Lake 

7.06 

5.96 

3.56 

5.31 

11 

+0.99 

Fish Lake 

16.33 

17.88 

16.45 

16.39 

33 

+2.34 

Big Lake 

13.24 

11.42 

12.00 

12.62 

25 

-1.89 

Goldwin Lake 

13.24 

13.91 

8.89 

11.07 

22 

+4.47 

Limes and Pond 

5.30 

2.48 

1.78 

7.08 

7 

-1.66 

Schelske Lake 

14.57 

12.91 

8.00 

11.29 

22 

+2.56 

Hust Lake 

38.85 

37.75 

31.11 

34.98 

69 

+4.35 

West Lake 

32.67 

30.79 

31.56 

32.12 

64 

-2.08 

Lawrence Lake 

5.74 

3.48 

5.78 

5.76 

11 

+3.59 

Eugene Pond 

7.06 

4.47 

4.89 

5.98 

12 

-2.37 

Woodworth Marsh 

4.86 

2.48 

1.78 

3.32 

7 

-1,31 


Totals 1034 +21.66 
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observations do not consistently follow this trend. Six of the water 
features observed by the SKYLAB scanner had computed areas greater 
than, six less than, and the remaining nine within the range of the 
respective LANDSAT observations. Since the SKYLAB observations 
occurred midway beOeen the two LANDSAT observations we have, for the 
sake of cotnpariso^f Assumed that the SKYLAB observed ponds should 
have had areas equal to the mean of the two LANDSAT observations 
(column d of Table 4). We then tabulated the departure of the actual 
SKYLAB observations from this mean and expressed the departure in terms 
of equivalent pixels (column e). An algebric summation indicated a 
cumulative deviation of only 21.7 pixels for the 21 ponds whose summed 
area was equivalent to a count of 1034 pixels. Thus it appeared that 
the SKYLAB observed individual pond areas varied randomly and that 
those ponds which were larger than expected tended to be balanced by 
those which were smaller than expected. He have examined the data 
including aerial photography and ground-truth photography in an attempt 
to explain these individual variations in the SKYLAB observations. We 
have not been able to relate these variations to differences in water 
quality, to the presence or absence of aquatic plants, to phenol ogical 
circumstances, or to differences in the near-infrared wavebands which 
were used to delineate surface water features*. We must conclude, 
therefore, that such factors were not contributory. Furthermore, it 
appears that for relatively small targets such as prairie ponds and 
lakes, the SKYLAB multi spectral scanner was not able to achieve as 
consistent a measure of area as was obtained with LANDSAT. We 
attribute this to the fact that the SKYLAB multispectral scanner 
utilized a conical scan. With such a system, problems of varying 
pixel overlap and data redundancy were inherent (Figure 23). The 
problems were further compounded when the data were converted into 
a straight line format using a nearest neighbor decision rule. The 
mapping of picture elements (pixels) of a conical scan-line into a 


*Water was delineated in the case of LANDSAT data with a 0.8 to 
1,1 urn waveband and, in the case of SKYLAB S-192 data, a 1.55 to 1.75 um 
waveband was used. 



FIGURE 23. RELATIONSHIP OF ADJACENT PIXELS IN THE GEOMETRY 
OF A CONICAL SCANNER. The figure depicts Increasing pixel 
overlap with displacement away from the ground track. 
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straightened scan-line must of necessity result in the discarding of a 
certain number of data pixels and slight dislocation of nearly all 
pixels from their true geographical position so that they can be fitted 
into a rectangular pixel grid. These actions have the net affect of 
slightly and systematically affecting areal representations and 
geometric fidelity. 

With regard to surface water area measurements* it is important 
to stress one other point. The pond sizes listed by the computer 
must in practice be termed “apparent size", because each pixel of 
data was examined and determined to be either totally water or not 
water. Many pixels lying on the perimeters of ponds and lakes undoubtedly 
contained some unrecognized and untabulated water. This caused the 
surface areas of virtually all water features to be underestimated. 
Percentage-wise, the errors were greater for the smaller ponds and 
for those of irregular shape {i.e., those having a high ratio of 
perimeter length to area). The very small ponds, of course, would 
not be recognized at all. In theory a pond must have filled or nearly 
filled the scanner's instantaneous-f1eld-of-view (IFOV) to be recog- 
nized*. The IFOV of the SKYLAB multi spectral scanner was 0.635 hectares 
in size. Recognition of a pond of this size would have been dependent 
upon whether the pond was wholly included in one digital sample or 
fractionally distributed over several samples. This would be governed 
by both the size and shape of the pond, by the frequency response of 
the scanner's electronics, and by the random occurrence of the pond 
with respect to the scanner's sampling grid (i.e. the occurrence of 
the pond with respect to a scan-line and the digital samples along 
that scan-line). In general, we feel that it is problematic whether 
ponds in the 0.6 to 2.5 hectare size class were recognized whereas 


*Since there was some variability in the radiance signals received 
from water targets (per the histograms of Figures 11 through 14), it 
is possible that a water target did not completely fill the scanner's 
instantaneous-field-of-view but was nevertheless dark enough to be 
classified as water. 


above 2.5 hectares it is most probable that they were recognized but 
not necessarily recognized at their full areal extent. The next chapter 
is devoted to the processing of multispectral data and particularly to 
the use of such data for detecting water elements smaller than the 
nominal resolution limits of the scanner. 


. ^ V 


CHAPTER 6 


PROCESSING OP MULTI SPECTRAL DATA FOR THE 
IMPROVED SPATIAL RESOLUTION OF WATER FEATURES 

As discussed in the previous chapter, signal level thresholding 
of a single near-infrared waveband of data is a simple and effective 
method for delineating surface water. With the advent of satellite 
programs, however, many workers in the field of remote sensing have 
been uncomfortable with the diminished resolving capabilities inherent 
in the operation of high altitude sensors, particularly scanners. Be- 
cause prairie ponds are frequently smaller than one-half hectare (1,2 
acre), it was apparent in both our LANDSAT and SKYLAB studies that many 
water features were not delineated by single channel thresholding. As 
part of these studies we have attempted to test a technique which takes 
advantage of the added information content of multiple spectral channels 
to estimate the proportion of materials present within a scanner's 
instantaneous-field-of-view (IFOV)*. The technique termed "proportion 
estimation" or "mixtures estimation" was first outlined by Horwitz et 
al. (1971) and further described by Nalepka et al. (1972). Before the 
LANDSAT and SKYLAB studies, the application of this technique was 
largely developmental in nature. Its use in these studies must be 
considered to be among the first attempts to test its applicability in 
a limited operational context. 

General Theory 

When the IFOV of a multi spectral scanner is large with respect 
to the scene objects being scanned, a single resolution cell may contain 


*The terms and "pixel" are often used interchangeably. 

However, the terms r, not synonymous for LANDSAT and SKYLAB data. 

The SKYLAB multi spectral scanner optics provide an IFOV of approxi- 
mately 79 X 79 meters while the data are sampled and digitized at a 
rate which equates a pixel to an area of approximately 72 x 72 meters. 
In the strictest sense, a description of the proportion estimation 
approach must make reference to the sensor's optical IFOV. In the 
actual processing and display of output data, however, reference will 
be made to the pixel or digitized sample- 
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a number of different material classes (i.e., the IFOV may be composed 
of a mixture of materials). The proportion estimation algorithm when 
applied to such data provides an estimate of the proportion of objects 
present within each of the scanner's IFOV's. 

A discussion of proportion estimation theory is beyond the scope 
of this paper, however, the essence of the technique can be described 
in geometric terms. Assume that a data set made up of two spectral 
channels, and contains three pure and unique materials -- 
A, B, and C. This situation can be depicted as in Figure 24 where 
the signature means for the three materials are shown in two-dimensional 
signal space. The signature simplex is the geometric figure formed 
by the lines connecting each pair of signature means. In the non- 
degenerate case, each pure signature is a distinct vertex of this 
simplex. If an unknown scene element (IFOV) consists of a mixture 
of all three materials, the signature of this material, X, lies within 
the simplex. An estimate of the pairwise proportion of pure materials 
constituting the unknown element is obtained by drawing a line from a 
vertex through the unknown signature to the opposite leg of the 
simplex. The inverse ratio into which each leg is divided defines 
the pairwise proportional composition of the unknown. In Figure 24, 
the unknown happens to lie at the centroid of the triangle, and its 
composition would be in the ratio of 1/3, 1/3, and 1/3 of materials 
A, B, and C, respectively. A case requiring special geometric inter- 
pretation is shown in Figure 25. In this instance the unknown, Z, 
lies outside or on the edge of the signature simplex. The unkown is 
determined to be made up of only materials A and C in the inverse 
ratio by which the simplex triangle's leg A-C is divided by a line 
drawn from Z orthogonally to that leg. If the unknown is quite 
distant from the signature simplex (described in terms of a distance) 
the algorithm is capable of designating the unknown as an alien object 
or an object composed of none of the simplex materials. 

Although the above description has been limited to three pure and 
unique materials in two-dimensional signal space, the concept is easily 
expanded to situations where many object materials exist in spectral 


FIGURE 24. GEOMETRIC INTERPRETATION OF MEANS OF 
SIGNATURE MIXTURES. In the case illustrated, the 
unknown, X, is a mixture of three pure materials 
— A, B, and G — which form the vertices of the 
signature simplex. 



FIGURE 25, GEOMETRIC INTERPRETATION OF ESTIMATE 
FOR A SPECIAL CASE. The unknown, Z, lying outside 
the signature simplex is a mixture of materials 

A and C. 
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hyperspace. In applyinq the alqorithm, however, it is necessary to 
observe two operational constraints. Firstly, at least n-1 spectral 
channels of information are required to satisfactorily estimate mixtures 
of n-object materials. Secondly, if the signatures for the materials 
in a mixture are similar or if one of them comes too close to a weighted 
average of the others, the estimates of the proportions may be poor. 

The latter condition is illustrated by Figure 26. Figure 26a shows a 
valid signal simplex for three signatures and two channels of data. 

Here covariance matrices interpretable in terms of loci of constant 
probability are shown. Figure 26b is a nearly degenerate signature 
simplex in which the vertex of one signature has come close to the 
weighted average of the other two signatures. A measure of what is 
"too close" is dependent upon the size and shape of the unit contour 
ellipsoid about the vertex or more specifically upon the covariance 
matrix. 

A Note About the Data Utilized For hulti spectral Processing 

In multispectral data processing, a necessary condition is that 
all channels of information utilized must be spatially registered. 

In accordance with the SKYLAB scanner design, all the even numbered 
SDO's as a group and all the odd numbered SDO's as another group should 
have been in registration within a group but one-half pixel out of 
register between groups. In evaluating several near-infrared wavebands 
as part of the single channel water mapping task (Chapter 5), it be- 
came apparent that misregistration existed not only between even and 
odd SDO groups but also within groups. This condition was discovered 
during an examination of several maps of two large lakes, each map 
generated by thresholding a different near-infrared waveband. If any 
two wavebands had beer spatially registered, all land/water interface 
pixels should have occupied the same geographic position in each of 
the respective threshold maps. The examination indicated, however, 
that for any two SDO's within either the even or odd numbered SDO 
groups, between 30 and 50 percent of a lake's peripheral pixels were 
randomly out of register by one or several pixels. 


(a) Signature Simplex with Unit 
Contour Ellipsoids 



FIGURE 26. GEOMETRIC CONFIGURATIONS FOR THREE 
SIGNATURES AND TWO CHANNELS 
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A similar lack of registration was independently discovered and 
more rigorously documented by Morgenstern et al., (1975) who used 
SKYLAB data collected over a Michigan test site. Band to band mis- 
registration could have occurred due to lags in the analog electronics 
or due to data displacements in the digital electronics either onboard 
the spacecraft or during ground processing. Braithwaite and Lambeck 
(1975) have shown that this source of misregistration was minor in most 
SDO’s. Morgenstern et al, (1975)s after examining data in both a conic 
scan format and a line straightened format, concluded that serious 
misregistration was created in the data by the scan-1 ine-straightening 
algorithm. 

The above stated registration errors posed potential problems to 
any multispectral processing and especially for proportion estimation 
processing. Consequently, we felt it inappropriate to use the scan- 
line straightened data which was currently on hand and which had been 
utilized for the single-channel water mappinn task (Chapter 5). lie, 
therefore, requested data for the same observation only in a conic 
scan-line format. These data were subsequently supplied to us by the 
Data Distribution Center of NASA/Johnson Space Flight Center. The 
remainder of this discussion will be devoted to the processing and 
analysis of this conic scan-line data. 

Proportion Estimation Processing 

In this phase of the study, the primary objective was to delineate 
open surface water in a mixutre of several scene materials. This 
should have made it possible to both improve the size estimates of 
larger ponds and to detect small ponds which would have been undetected. 
Since models currently used for estimating waterfowl production 
utilize pond numbers, we have emphasized the detection and enumeration 
of ponds rather than their areal measurement. Secondarily, we were 
interested in detecting wetland components which were characterized by 
marsh vegetation, the canopies of v/hich largely occlude standing 
surface water. Such marsh conditions are usually peripheral to open 
water and are often too small to be delineated in whole pixel recoonition. 
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Before processing began, the multi spectral data {in a conic scan 
format) was evaluated by visually assessing electronic screening imagery 
produced from the data tape. Table 5 presents a summary of that assess- 
ment while Figure 27 presents examples of both good and poor quality 
imagery. Note that in terms of target contrast, only those bands in 
the near or reflective infrared were of good quality. All visible 
and thermal -infrared bands were of lesser quality in terms of contrast 
and other forms of electronic noise. This condition can perhaps be 
explained by the time of day during which the observation occurred — 
12:56 GMT or 06:19 local solar time. This places the observation at 
about two hours after sunrise, at a time when the solar altitude was 
less than 20°. It is conjectured that the near- infrared detectors 
produced a cleaner output signal under these conditions than did the 
visible detectors because; (1) with the longer atmospheric path 
occasioned by the low sun angle, infrared radiation was scattered 
considerably less than was visible radiation, (2) the bandwidths of 
the infrared detectors were broader than were the bandwidths of the 
visible detectors, and (3) a vegetation dominated terrain has generally 
a higher reflectance in the infrared than in the visible. Furthermore, 
the SKYLAB multispectral scanner was designed for optimal performance 
at solar altitude angles in excess of 30° (National Aeronautics and 
Space Administration 1973). The poor contrast in the thermal infrared 
is again attributed at least in part to the time of the observation. 
Generally, during a brief period shortly after dawn and again after 
sunset, temperature differences between most terrain objects will be 
muted. This is due to the warming effect of the sun and the differing 
heat capacities of terrain features. 

Figure 28 further illustrates characteristics of the data utilized. 
In this case, a one-percent systematic sample was extracted from the 
data set. The figure indicates the data value range of 95 percent of 
the sampled pixels. Generally the samples had an approximate Gaussian 
(normal) distribution between the extremes shown. For the observation 
of 12 June 1973, the data values were generally depressed and lacked 
dynamic range. These depressed data values can largely be attributed 


TABLE 5. EVALUATION OE SKYLAB MULTISPECTRAL SCANNER ELECTRONIC SCREENING 
IMAGERY EOR DATA OF 12 JUNE 1973 COLLECTED OVER EASTERN NORTH DAKOTA 


BAND A range Cum) 

SDO 

IMAGERY CONTRAST 

LOW FREO. NOISE 

NOISE BANDING* 

1/f NOISE 

COMMENTS 

1 0.41-0.46 

22 

very 

poor 

yes 

yes 

no 


2 0.46-0.51 

18 

poor 


yes 

yes 

no 


3 r 0.52-0.56 
~Lo.52-0.56 

1 

poor 

to fair 

yes 

yes 

no 

some sync, dropout 

2 

poor 

to fair 

yes 

yes 

no 

some sync, dropout 

4 ro. 56-0. 61 
LO. 56-0. 61 

3 

poor 


yes 

yes 

no 

some sync, dropout 

4 

poor 


yes 

yes 

no 

some sync, dropout 

5 ro. 62-0. 67 
“L0.62-0.67 

5 

very 

poor 

yes 

yes 

no 


6 

very 

poor 

yes 

yes 

no 


ro. 68-0. 76 

^ “Lo. 68-0. 76 

7 

good 


yes 

no 

yes 


8 

good 


yes 

no 

yes 


. ro. 78-0. 88 
LO. 78-0. 88 

9 

good 


no 

no 

no 


10 

good 


no 

no 

no 


8 0.98-1.08 

19 

good 


no 

no 

no 


9 1.09-1.19 

20 

good 


no 

no 

no 


10 1.20-1.30 

17 

good 


no 

no 

no 


n _r 1.55-1.75 
LI. 55-1. 75 

11 

good 


no 

no 

no 


12 

good 


no 

no 

no 


r 2.10-2.35 
L2. 10-2.35 

13 

good 


yes 

no 

no 

some sync, dropout 

14 

good 


yes 

no 

no 

some sync, dropout 

in iS~ 10.20-12.50 
^ “LlO. 20-12. 50 

15 

ve ry 

poor 

yes 

no 

no 

some sync, dropout 

16 

very 

poor 

yes 

no 

no 

some sync, dropout 

13-1 10.20-12.50 

21 

very 

poor 

yes 

yes 

no 



*Noise banding is noise occurring In phase with the scan frequency. In these particular cases it was 
evidenced by two cycles of alternate dark and light bands which occurred throughout the imagery 
parallel to Che ground crack. 
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0.46-0,51 0,56-0.61 0.68-0.76 0.98-1.08 1.20-1.30 2.10-2.35 

WAVEBAND (miii) 


FIGURE 28. DYNAMIC RANGES OF SKYLAB MULTISPECTRAL SCANNER DATA COLLECTED OVER EASTERN 
NORTH DAKOTA ON 12 JUNE 1973 AND, FOR COMPARISON, OVER SOUTHERN MICHIGAN ON 5 AUG. 1973. 
The ranges are based on a 1% systematic sample taken throughout the data. 
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to the low sun angle and the generally poor illiunination conditions 
which existed. For ccjnparison, the figure also includes similar data 
samples obtained from a SKYLAB scanner observation made over southern 
Michigan on 5 August 1973 at 15:00 GMT or about 09:30 local solar time. 
Note that these latter samples have much broader dynamic ranges in 
nearly all spectral bands and that the ranges are not depressed. 

The actual implementation of proportion estimation processing 
involved as a first step the securing of spectral signatures for object 
materials occurring in the observation scene. Multispectral signatures 
extracted from actual scene elements (training sets) for the SKYLAB 
observation of 12 June 1973 are shown in Figure 29. Only even numbered 
SDO's were utilized and, in addition, data in the 0.41- to 0,46-um 
and 10.20- to 12.50-um wavebands were discarded because of poor 
quality. In selecting the training sets, care was taken to pick 
resolution elements that were pure in their constituency. In order 
to obtain representati ve samples, however, the signatures were obtained 
by combining several training sets which consisted of like materials. 

For example, the water signature represents a combination of several 
ponds and lakes which ranged in water quality from the relatively 
clear to the moderately turbid. As a result, the sianature for water 
in this instance has a larger standard deviation than would normally 
be expected. Similarly, other signatures represent a variety of field 
and marsh situations. Note from the figure that the signatures are 
not well differentiated in the visible wavebands -- a further mani- 
festation of poor contrast as observed for these bands in the electronic 
screening imagery. 

The deep marsh signature was obtained from several communities 
of bulrushes ( Scripus spp.) and cattails ( Typha spp.). Bulrushes 
commonly occur in solid stands and frequently in association with 
cattails. Cattails occur less frequently in solid stands and conse- 
quently were not as dominant in the deep-marsh composite sionature 
as were the bulrushes. Signatures for the shallow marsh class were 
obtained from plant associations of whitetop ( Scolochloa festucacea ) 
and sedges (Carex, spp.). Whitetop is a tall, lush marsh grass that 
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FIGURE 29. SPECTRAL SIGNATURES USED IN THE PROCESSING OF SKYLAB MULTISPECTRAL SCANNER 
DATA COLLECTED 12 JUNE 1973 OVER EASTERN NORTH DAKOTA. The signature mean and plus and 

minus one standard deviation are shown. 
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often grows 1n solid stands and, in drier years, is cut for hay. 

Haying stimulates the growth of whitetop allowing the species, in 
the long term, to become dominant over the sedges. On the other 
hand, pasturing tends to suppress the whitetop and allows sedges to 
become dominant. Both of these marsh classes frequently occur in 
narrow concentric rings surrounding central areas of open water. 
Continuous, large expanses of these communities are less common. It 
is this latter situation which must be selected as a training set and, 
as a result, the selection was often limited. 

The small grain signature represented those grains typically grown 
in this area of Morth Dakota — hard red wheat, durum wheat, barley, 
and oats. These grains, all springplanted, were in an early growth 
stage during the observation on 12 June. Grain fields throughout 
the region were represented by a wide range of phenol ogi cal stages 
because of variations in planting dates. Some fields had only recently 
been planted, thus bare soil was predominant. Other fields had been 
planted for up to one month and had reached a stage where a nearly 
closed, homogeneous canopy existed. It is this latter condition 
that is represented in the small grain signature. The "idle" signature 
represents primarily grasslands (i.e. , pasture, open range, and/or 
native prairie situations). In general, grasslands in early June 
were represented by a predominance of above ground, dead, herbaceous 
standing or matted biomass remaining from the previous year's growth of 
grasses and forbs. Few green plants had as yet emerged from this 
standing or matted dead biomass. 

The signa^■ure5 selected v;ere evaluated with an automated 
statistical analysis program. The purpose of the analysis was to 
determine whether the position of the signatures in mul ti spectral 
hyperspace permitted a meaningful mixture estimation. If, for example, 
three signatui^e means. A, B, and C, were in a line, then one, would 
have no way of knowing whether a data point between B and C was a 
mixture of B and C, a mixture of A and C, or a mixture of all three 
(an approximation of this situation was illustrated previously in Finure 
26b). The results of our analysis of a set of six, five, and four 
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signatures respectively are presented in Table 6. In the first 
case, where the full complement of six signatures was tested, the 
separability of four of the six signatures was poor and for water 
and bare soil was only moderate. In the second iteration the shallow 
marsh signature was eliminated, but water's separability did not im- 
prove. Finally the set was pared to four signatures, at which time the 
separability of all signatures, but particularly of water, improved 
drastically. Our rationale in discarding the deep marsh signature 
was that this signature probably contained some water along with 
vegetation particularly due to the timing of the observation when 
water levels were high and the new year's vegetative growth had not 
■Pully developed. In eliminating the shallow marsh signature we felt 
that its characteristics were duplicated by either the small grain 
signature (if a flush of new growth had occurred) or by the "idle" 
signature (if the present year's growth was still masked by the previous 
year's dead biomass). 

As a result of the analysis described above, we chose to perform 
proportion estimation processing using four signatures - water, bare 
soil, small grain, and idle. For this set of signatures the small grain 
signature had in essence become a surrogate for all green herbaceous 
vegetation. {Woody plants at any significant scale were not present.) 

The fact that the estimation was to be done without the deep marsh 
signature also meant that some of these marsh components would be 
recognized and tabulated as open water. We felt it better to bias 
the error in this direction as compared to drastically underestimating 
water, particularly pond numbers. In fact, the transition between 
open water and closed stands of emergent plants is often a continuum 
which can vary both with the season and with position within a wetland. 

As a result, the delineation between open water and an emergent marsh 
community is often a judgment decision. 

The proportion estimation algorithm was applied to multispectral 
information in a conical scan-line format. The output of this processinn, 
for purposes of display and statistical analysis, was then scaled and 
converted into a scan-line straightened format. The line straigr.iienino 


TABLE 6. ANALYSIS OF SIGNATURE SEPARABILITY - A Test for the Uniqueness 
of a Signature for Use in Proportion Estimation Processing 


Number of 
Signatures Used 
for Processing 

Water 

Deep 

Marsh 

Distan 

Shallow 

Marsh 

ce* 

Bare 

Soil 

Small 

Grain 

Idle 

6 Signature Set 

2.2788 

0.7810 

0.4534 

2.4743 

1.2150 

0.9668 

5 Signature Set 

2.3477 

1.6264 


2.8670 

4,6227 

1.4359 

4 Signature Set 

7.3618 



3.2506 

4.8842 

2.7381 


* The distance in standard deviation units of the signature mean 
from the weighted average of the remaining signatures. 
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algorithm was similar to that which had been used by the Earth Resources 
Production Processing System of NASA/ Johnson Space Center except that 
compensation for the earth's rotation was not performed. 

Results 

For the proportion estimation computation, the output was a set 

of proportions for each pixel. A water recognition map generated from 

this output is shown in Figure 30. For comparison a map generated with 

the single-waveband thresholding algorithm is shown in Figure 31. The 

data in both figures have been scaled (so that the digital pixels have 

a height/width ratio compatible to a line printer) and converted to a 

2 

line straightened format. Both maps represent the identical 286 km 
(110 square mile) area. In the proportion estimation map, the symbol 
density is related to the proportion of water estimated for that 
pixel. In order for the map to accurately portray the scene, certain 
percentage or acceptance limits were determined for the output of the 
algorithm. For example, it seemed appropriate to count pixel values 
of 0.81 and above as totally water. This procedure tended to account 
for the likelihood that a value close to a signature mean (i.e., close 
in terms of the probability contours) may in fact have been a pure 
sample related to that mean. Similarly, pixels showing less than 0.40 
water were assumed to be false alarms (i.e., nonwater pixels classified 
as water), and they were excluded from any consideration as surface 
water. These limits were established after an examination of a small 
portion of the processed data and a comparison with mul tispectral 
scanner data and photography collected by supporting aircraft. 

In general, a detailed comparison of the classification maps and 
related imagery indicated that proportion estimation significantly 
improved both pond shape definition and the recognition of smaller 
water features which otherwise would not have been detected. 

In the proportion estimation processing of LANDSAT data, a 
difficulty encountered was the inability to adequately delineate 
alkaline lakes. Such lakes are scattered throughout prairie areas 
of non or poorly integrated drainage but particularly in glacial 




Fish Lake 


FIGURE 30. WATER RECOGNITION OBTAINED BY THE USE OF THE PROPORTION 
ESTIH.\TION ALGORITHM APPLIED TO SKYLAB SCANNER DATA COLLECTED 12 
JUNE 1973. The symbol density is related to the proportion of water 

estimated for that pixel. 
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FIGURE 31. WATER RECOGNITION OBTAINED BY THE USE OF THE SINGLE 
CHANNEL THRESHOLDING ALGORITHM APPLIED TO SKYLAB SCANNER DATA 
COLLECTED 12 JUNE 1973. The decision criteria is such that each 
pixel has been classified as either totally water or not water. 
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outwash areas of the Missouri Coteau. A proportion estimation water 
map which closely matches the areal coverage of our SKYLAB map (Figure 
30) and which was produced from LANDSAT data collected 25 days after 
the SKYLAB observation is included here as Figure 32, Both Alkali 
Lake and Norden Lake in the LANDSAT map were only partially recognized. 
Both lakes at the time contained large amounts of suspended solids 
and possibly precipitated alkali bottom sediments. These conditions 
are especially prevalent during periods of low water as was the 
situation during the summer of 1973. The anomalous detection of these 
lakes, we had felt, could have been improved had there been additional 
spectral information available beyond what was available with the 
LANDSAT multi spectral scanner. In particular, we had thought that 
the additional near-infrared wavebands and especially the 1.55- to 
1.75-|jm waveband available on SKYLAB would provide improvements in 
the recognition of these anonialous lake features. A comparison of 
Figure 30 and 32 indicates that a dramatic improvement was, in fact, 
realized. We are unable, however, to determine in a quantitative 
sense how accurately the areas of the two lakes were estimated since 
no current low altitude planimetric data for the lakes were available. 

Quantitative comparisons were made of certain other surface 
water features which coincided with an aircraft data transect located 
in the vicinity of the Woodworth Station. The strip maps of Figure 33 
were plotted from processed SKYLAB scanner data while still in a conic 
scan-line format. The data were quantitatively analyzed in this format 
in order to avoid errors which would be introduced by the scan-line 
straightening algorithm. The map symbols used are plotted at the 
areal centroid of the respective pixel, the symbol size being equivalent 
to the proportion of water found in that pixel. Note that the symbols 
are not equally spaced. This is due to the overlapping of adjacent 
pixels in a direction orthogonal to the ground track of the satellite 
and in direct proportion to the distance of the pixel from the ground 
track (see Figure 23). 

Figure 33a resulted from the thresholding of a single waveband of 
data (1.55 to 1.75 urn). Figure 33b is a proportion estimation map 
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(b) Water Recognition by Proportion Estimation Processing of SKYLAB 
Multispectral Scanner Data Collected 12 June 1973. The minimum frac- 
tional acceptance limit for water was 40% per pixel. 



(c) Water Recognition by Proportion Estimation Processing of SKYLAB 
Multispectral Scanner Data Collected 12 Tune 1973. The minimum frac- 
tional acceptance limit for water was 31% per pixel. Pixel water 
fractions in the range 31 thru 39% are Indicated by small open symbols. 

FIGURE 33. WATER RECOGNITION IN THE VICINITY OF WOODWORTH STATION, 
NORTH DAKOTA. The proportion estimation recognition maps (b and c) 
differ only in the minimum fractional acceptance limit used in plot- 
ting the map. In both of the proportion estimation maps, symbol size 
is related to the percent of water detected in that pixel. The symbols 
labelled "E" are confirmed commission errors. Scale; 1/62,500. 
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vvhich uses a minimum acceptance level of 40 percent ('•’.e., pixels 
tabulated as containing less than 40 percent v/ater v^ere assumed to 
contain no vjater at all). In this map, it is appart;i;t that a major 
improvement in both the detection of small ponds and the detection 
of peripheral pond features was realized. The map contains two con- 
firmed commission errors at locations "E". Many other smaller ponds 
were not identified as may be attested by a comparison v/ith the 
electronic aircraft imagery of Figure 34. 

Referring again to Figure 33b, the smallest confirmed recognized 
water feature was Sargent Pond. Figure 33c is also a proportion map 
similar to the previous map except that an acceptance level of 31 
percent was utilized. (Pixels in the 31- to 39-percent range are 
indicated by small open symbols. The map differs only from the 
previous map with the inclusion of these symbols.) Clearly many of 
these added symbols do represent vrater, but it is equally evident 
that many also v/ere com;„ission errors as at locations labelled "E". 

It was because of the proliferation of these errors that the minimum 
acceptance level was set at 40 percent. In addition to these maps 
which can only qualitatively illustrate a range of proportions, the 
areal extent of each of the several ponds and lakes named in Figure 
32b, was calculated based upon the exact proportion of water listed 
in the algorithm output stream (and excluding pixels containing ss 
than an estimated 40 percent water). 

Before describing the results, ho'wever, vje present here a brief 
discussion of the areas which vjere assigned to each pixel. The IFOV 
of the SKYLAB multi spectral scanner was 80 x 80 m but, since the 
scene was overscanned in the direction of the satellite's velocity 
vector and since the data v/ere oversanipled in a direction orthogonal 
to this vector, there v/ar- overlap in the ground patch covered by 
adjacent pixel samples. In calculating the v/ater area of a pond or 
lake, one needs to consider the actual area viev/ed by each pixel. In 
other v/ords, if a pond smaller than SO x 80 in is contained v/ithin 
one pixel, the pond area is 50 percent of 80 x 80 m (i.e., the IFOV) 
and not 50 percent of the smaller effective area. Now if this same 


vjT'tiOnUCIBILITy OF THE 
• FAGE is P(X)R 



(a)j.O- to 1.4- viin Electronic Imagery — 12 May 1973 from 1370 m Altitude. 
(The imagery appears skewed due to a pronounced aircraft crab angle.) 


(b)l.O- to 1.4- ym Electronic Imagery — 12 Aug. 1973 from 1370 m Altitude. 


FIGURE 34. AIRCRAFT MULTISPECTRAL SCANNER VIDEO COLLECTED ON A TRANSECT OVER WOODWORTH 
STATION, NORTH DAKOTA. Approximate along track scale: 1/26,000. 
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pond was seen in the overlap area of two adjacent pixels it would 
be inaccurate to usc the 80 x 80 m area for each pixel since some 
portion of the pond would be counted twice. To account for problems 
of this sort, allowances for adjacent water pixels v/ere made when 
calculating estimated pond and lake areas. (Any pixel containing 
water may have been overlapped by up to a maximum of six other pixels 
containing water - as many as two in the direction of the satellite's 
velocity vector and as many as four orthogonal to this vector.) 

In addition to the pond and lake areal measurements obtained with 
SKYLAB scanner data, more precise measurements of area were also 
obtained with the single channel thresholding algorithm applied to 
low altitude aircraft scanner data. These latter data were used to 
compare the SKYLAB proportion estimation results against. In Table 
7 the results are presented for comparison. Note that neither of the 
two aircraft observations coincided with the satellite observation 
but preceded and followed it by 31 and 60 days respectively. From the 
comparisons, six of the SKYLAB observed water features were reasonably 
close in area to the aircraft observations (Big and Fish Lakes, Pond 
8-14, Glen's Pond, Sargent Pond, and Koening Pond). One water body 
was underestimated in area (Clark Lake), The remaining six water 
features v/ere overestimated in area. Thus it would appear that the 
proportion estimation results tended to either approximate or somewhat 
overestimate the actual areas of ponds and lakes. This would seem 
to bear out our earlier suspicions that areas of emergent marsh 
plants, may in certain cases, be tabulated as open water areas. In 
particular, Goldwin Lake, Carl's Pond, and Woodworth Marsh are known 
to have had extensive peripheral deep marsh areas, and the areas of 
all of these features were notably overestimated. We feel that efforts 
to effectively estimate open surface water and more especially areas 
of marsh communities were compromised by inadequate data quality in the 
visible bandwidths. As noted previously, this lack of quality appears 
to have been occasioned by the early hour of the SKYLAB observation. 

Finally, we present the results of areal water tabulations for 
2 

the 286 km scene of Figures 30 and 31. These tabulations were 


TABLE 7. TABULATION OF AREAL MEASUR^ffiNTS OF OBSERVED POND AND LAKE FEATURES 


Computed Areas (Hectares) 

Aircraft Observation SKYLAB Observation Aircraft Observation 

of 12 May 1973 (1) of 12 June 1973(2) of 12 Aug. 1973 (1) 


Liraesand Pond 

Not Available 

6.15 

4.49 

Goldx-fin Lake 

16.13 

23.59 

13.78 

Pond 9-1 

2.41 

3.20 

2.67 

Big Lake 

17.56 

18.90 

16.81 

Fish Lake 

24.42 

24.76 

23,53 

Clark Lake 

11.43 

8.36 

9.04 

Pond 8-lA 

00 

o 

1.15 

Dry 

Glen's Pond 

2.43 

2.14 

Dry 

Sargent Pond 

0.34 

0.28 

0.06 

Sadie Pond 

2.36 

2.90 

2.53 

Koenig Pond 

1.25 

0.80 

Dry 

Carl's Pond 

1.10 

3.43 

Dry 

Woodworth Marsh 

5.75 

6.88 

3.06 


Notes: (1) Pond areas were computed by thresholding a 1.5- to 1.8- pm waveband of 

scanner data. Observation was made from an altitude of 1370 m. 


(2) Pond areas were computed by proportion estimation processing of SKYLAB 
multispectral scanner data. 
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accomplished using the computer program described in Appendix A. 

This program is not able to take account of differences in area 

represented by the scanner's IFOV and the digital sample (pixel) 

nor of pixel overlap. Therefore, the results can be compared only 

in a qualitative sense. Table 8 presents a comparison of numerical 

2 

and area tabulations of water features for the 286 km test area 
using both the threshold and proportion estimation algorithms. 

This comparison indicates that the total number of ponds and lakes 
delineated by proportion estimation processing was 189 percent of 
the total number obtained with the single channel threshold algorithm. 
The results of proportion estimation processing of LANDSAT data for 
the identical scene are also tabulated in Table 8. It is significant 
to note that proportion estimation processing of both SKYLAB and 
LANDSAT data delineated nearly an equal number of water features. 

The distribution of pond numbers within size classes particularly the 
smaller classes differs because the basic pixel sizes were different*;' 
this difference allowed the ponds to be clustered somewhat differently 
into the size classes shown in the table. 


*Pixel size for the LANDSAT data was 57 x 79 m. The SKYLAB pixel 
after we had line straightened and scaled the data was 58 x 71 m. 



TABLE 8. COMPARISON OF TABULATIONS OF PONDS AND LAKES 


SKYLAB Observation of 

12 June 1973 


LANDSAT 

Observation of 7 June 1973 

Computed Using the 
Thresholding Algorithm 

Computed Using the 
Proportion Estimation Algorithm 

Computed Using the 
Proportion Estimation Algorithm 

Area (Acres)* Frequency 

Area (Acres)* 

Frequency 

Area (Acres)* 

Frequency 

0.25 to 0.50 


0.25 to 0.50 

109 

0.25 to 0.50 

72 

0.51 to 1.00 

56 

0.51 to 1.00 

35 

0.51 to 1.00 

71 

1.01 to 2.00 

26 

1.01 to 2.00 

52 

1.01 to 2.00 

62 

2.01 to 3.00 

11 

2.01 to 3.00 

33 

2.01 to 3.00 

31 

3.01 to 4.00 

11 

3,01 to 4.00 

18 

3.01 to 4.00 

19 

4,01 to 6,00 

16 

4.01 to 6.00 

24 

4,01 to 6.00 

21 

6.01 to 8.00 

10 

6.01 to 8.00 

21 

6,01 to 8.00 

20 

8.01 to 10.00 

9 

8.01 to 10.00 

6 

8.01 to 10.00 

11 

10.01 to 15.00 

10 

10.01 to 15.00 

15 

10.01 to 15.00 

15 

-.5.01 to 20,00 

13 

15.01 to 20.00 

10 

15.01 to 20.00 

7 

20.01 to 25.00 

5 

20.01 to 25.00 

7 

20.01 to 25.00 

5 

25.01 to 30.00 

0 

25.01 to 30.00 

5 

25.01 to 30.00 

3 

30.01 to 40.00 

3 

30.01 to 40.00 

2 

30,01 to 40,00 

4 

40.01 to 50.00 

6 

40.01 to 50.00 

8 

40.01 to 50.00 

6 

OVER 50.00 

U 

OVER 50.00 

19 

OVER 50.00 

14 

TOTAL 

193 


TOTAL 364 

TOTAL 

361 


^Although Metric un’" cs of measure are general ly used throughout this text, the computer software In current 
use clustered water features according to English units of areal measure (acres). 
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APPENDIX A 


DESCRIPTION OF COMPUTER PROGRAM 
FOR GENERATING WATER BODY AREA AND PERIMETER STATISTICS 

The software program, APSTAT (Area, Perimeter Statistics), is 
designed to delineate bodies of open surface water and to generate 
statistics (area, perimeter, and shape factor) on these water bodies, 
APSTAT has evolved from an older program, MAPLKS, also developed by 
the Environmental Research Institute of Michigan. Both programs 
utilize the same decision criteria for classifying a grouping of 
pixels as a pond or lake. Both programs compute the area of any 
specific closed feature (i.e., a pond or lake) one line at a time, 
summing the areas in each scan line to determine the area of the 
specific feature. MAPLKS takes account of increases in the spatial 
field of view (and thus an increase in pixel area) as scan angles 
increase from the nadir position. APSTAT does not take account of 
variations in pixel area but instead assumes that all scene pixels 
are of a constant area. Thus APSTAT is suitable for use with 
satellite data where scan angles subtend no more than 12 or 15 
angular degrees. Specifically, APSTAT has been designed to operate 
on data collected using LANDSAT and S'<''LAB/EREP scanners both of which 
have approximate lateral fields of view of 11”. 

Although APSTAT was developed primarily for the tabulation and 
analysis of pond and lake occurrence, application of the program 
need not be limited to water bodies. Hs statistical enumerating 
capability can be applied to any scene feature having a closed peri- 
meter such as forests or agricultural crops provided there exists 
digitized data in which the scene feature of interest can be de- 
lineated on the basis of a discrete voltage range in a single data 
channel (as, for example, a classification tape or water features 
which exhibit uniquely low radiance values in a near infrared waveband). 

Specifically, APSTAT can identify any class of data which may be 
defined by one of two modes of operation defined by the following level 
slicing algorithms: 
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1. 0 £ (pixel voltage count on channel ICODE) < VHI» 

the pixel is assigned to the scene class of interest. 
Otherwise: the pixel is assumed not to belong to that class. 

2. If: V0LTS(([i) ^ (pixel voltage on channel ICOD) ^^VOLTS(l), 

the pixel contains PC(<|))% of the class of interest. 

If; VOLTS(l) _< (pixel voltage on channel ICODE) £V0LTS(2), 
the pixel contains PC(1)% of the class of interest. 

If: V0LTS(2) ;< (pixel voltage on channel ICODE) ^V0LTS(3), 

the pixel contains PC(2)% of the class of interest. 

If: V0LTS(3) (pixel voltage on channel ICODE) ^V0LTS(4), 

the pixel contains PC(3)% of the class of interest. 
Otherwise: the pixel does not belong to that class. 

Where ICODE, VOLTS(O). . ,V0LTS(4) and PC(1), . .PC(3) are 
user specified. 

The first algorithm is the normal (default) mode of classification 
The second is the "proportion estimation" mode and must be called by 
the user. The proportion estimation mode assumes the availability of 
multiple channel tape with each material present in the scene re- 
presented by one channel — the integer scale on the channel being 
representative of the proportion of that material present within each 
pixel. Only one channel (i.e., one material class) of a proportion 
estimation tape can be processed at a time. 

APETAT also incorporates the following features; 

(a) The program will accommodate up to 1040 points per linear 
scan line. 

(b) The program can accommodate multiple channel input (to a 
maximum of 13) but must operate on only one of these channels 

(c) A maximum of 3250 lines of data may be processed each time 
the program is run. 

(d) The program lists the position of each identified feature 
by the last (highest) scan line on which it appears and the 
last (highest) point of the feature on that scan line. 
Optionally, the feature may also be identified by the 
latitude and lonoitude of that point. 


Algorithm Descriptions 
1 ) AREA 

In the normal (default) mode of area calculation, the area of 
a feature (such as a lake) is equal to the number of feature pixels 
times the area of a pixel (always constant). In mathematical terms: 

Feature Area = Number of Feature Pixels x Pixel Length x Pixel Width 

where pixel length and width are user specified. 

In the proportion estimation mode 


Feature Area = Pixel Length x Pixel Width 
3 

X (number of feature pixels which are PC(i)%) 
i=0 

2) PERIMETER 

The following examples illustrate the definition of perimeter: 



Each pixel identified as the specified feature is shown as an X. The 
calculated perimeters are shown by the solid dark lines, and the 
arrow points to line segment whose lengths are counted twice. The 
following characteristics are noted: 

(a) Perimeter measurements are made from pixel centers, not from 
pixel margins, 

(b) Perimeter calculations cannot be performed in the proportion 
estimation mode. 
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ORIGINAL PAGE IS POOR 


(c) The perimeter of features which overlap a scene (or 
universe of interest) boundary will include a false 
boundary. For example: 



Although the above lake continues beyond the scene 
boundary, the perimeter is calculated as the sum of the 
line segments shown. 

3) SHAPE 

The shape factor is defined as the ratio of the perimeter to 
the square root of the area normalized to one {1.0) for a round lake. 
The larger the shape f.;;rLor, the more irregular is the lake's shore** 
line. Mathematically the shape factor may be represented as follows: 

Shape Factor = £er,<jeter ^ j_ 

/Area 2/ii~ 

where is the normalization factor. 

2)'^TT 

Because perimeter measurements are made from pixel centers, not 
from pixel margins, shape factors for lakes of size less than 10 or 12 
pixels may not be valid. 
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4) LATITUDE/LONGITUDE 

The latitude and longitude of a feature are defined c^s the 
latitude and longitude of the last (highest numbered) point of that 
feature on the last (highest numbered) scan line on which the 
feature occurs. Transformation coefficients to convert from line 
and point coordinates to geographic coordinates are input by the user 
and are obtained by a linear regression analysis external to tnis 
program. (Because the coefficients currently in use are first order 
terms and because lines of latitude converge toward the poles, the 
areal extent over which one set of coefficients can be applied 
should be no greater than 25 x 25 nautical miles. In order to 
preclude this problem in future usage, the use of the Universal 
Transverse Mercator coordinate system is contemplated. 

5) CLASSIFICATION 

Recognized features (lakes) are delineated by certain rules of 
classification algorithms. The following examples illustrate the 
rules by which pixels are grouped into an areally limited feature 
(lake): 
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Each letter (Q, R, S, T, Vj W, )(, Y» Z) represents a pixel. Groups 
of the same letter represent regions which would be classified to- 
gether as forming one lake* (if this is the feature being recognized). 
Briefly the rules used to group the pixels are as follows: 

(a) A continuous series of pixels on a scan Tine will be 
considered as constituting a lake as in lake Z, or a part . 
of a lake as in Y, V. W, and S. A discontinuity in water 
pixels may be bridged and counted as water only if the 
bridging option is utilized (see 6 below). 

(b) Pixels in a subsequent scan line are again grouped if (a) 
above applies, and the water segment will be 1 inked to 
pixels in the previous scan line if any* some, or all of 
the pixels are vertical Ty or diagonally adjacent to any, 
some or all of the water pixels in the previous scan line. 
Lake R illustrates the simple case of a lake consisting of 
one pixel in each of two adjacent lines, the pixels being 
vertically adjacent. Lake Q is a two pixel lake consisting 
of one pixel in each of two adjacent lines, the pixels being 
diagonally adjacent. By means of this algorithm it is 
possible for several arms of one lake to be connected as 

in Lake X. 

6) BRIDGING 

Pixels not original ly identified as water by level slicing may 
be redefined as water by the bridging algorithm, where the user 
specifies; 

(a) the maximum number of adjacent non-water pixels on a single 
scan line which may be bridged and 

(b) the minimum number of adjacent Abater pixels on a single 
scan line which must be in the water intervals on each side 
of any non-water segment for the non -water segment to be 
bridged. 

NOTES; (i) Bridging fs not possible when the proportion 

estimation mode of classification is specified. 
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(ii) Bridging only occurs on one scan line, never between 
scan lines. 

7) AREAL STRATIFICATION 

The user my wish to stratify the statistical results. The 
program is capable of sequentially processing any number of strata 
(scenes). However, a scene must consist: of no more than ten 
rectangles each of which are specified in terms of beginning and 
ending scan lines and beginning and ending points. For example: 


^ POINTS — ^ 



STRATUM 


“a" 


In the above example stratum '‘a" consists of one rectangle, 
stratum "b“ of an aggregation of four rectangles, and stratum "c“ of 
two rectangles. The size of each stratum is limited only in the 
number of rectangles which may be grouped to comprise a stratum (ten), 
and by the maximum number of scan lines (3250) and maximum number of 
points within a line (1040) which the program can handle. Appendix B 
exemplifies a typical APSTAT statistical listing for two strata of data 




, : APPENDIX. B 

TABULATION OF POND AND LAKE STATISTICS 

This appendix incTudes a tabulation of recognized ponds and 
lakes identified throughout the study area using SKYLAB/EREP multi - 
spectral scanner data collected on 12 June 1973. Based on physio- 
graphic differences, the study area was divided into two strata— . 
the coteau stratum and the drift plain stratum respectively. The 
drift plain stratum was further subdivided into units labeled 
"driTt plain west" and "drift plain east". This subdivision was 
necessary because the raw data were divided between two separate 
computer compatible tapes which could not be abutted for pro- 
cessing. The following diagram and table indicates the vertices 
in geographic coordinates of the several strata and substrata. 




STUDY AREA AND STRATA GEOGRAPHIC COORDINATES 


Vertex Code 


A 

B 

C 

D 

E. 

F 

G 

H 

I 

J 

K 

L 

M 

N 

0 

P 

Q 

R 


Vertex Coordinates {degrees) 


Latitude 


47.595 

47.356 

46.972 

47.210 

47.069 

47^453 

47.129 

47.147 

47.201 

47.224 

47.265 

47.288 

47.317 

47.341 

47.350 

47.374 

47.396 

47.427 


Longitude 


99.401 

98.431 

98.648 

99.618 

99,043 

98.826 

99.009 

99.080 

99.050 

99.144 

99.121 

99.215 

99.199 

99.293 

99.287 

99.382 

99.369 

99.495 
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compoter-tabulatT^^ list each, water feature 

recognized in the stratum by enumerating the location and size of 
the feature i The location is defined using each of two coordinate 
systems based (1) on a scan Tine and point number scheme and (2) 
on a more conventi orial gedgraphi c or 1 ati tude and 1 ongi tude scheme . 
The convention used to reference each water body has been to locate 
the body at the position of the last (highest numbered) scan line 
with at least one pixel in the water body and the last (highest 
,■'■;■■■■■ ■ numbered) water body pixel of that scan line. The water feature. area 
■■■ is. given in terms of both English (acre) and metric (hectare) units. 
Following this enumeration , the assigned name of the scene i s listed 
along with a description of the scan lines and points which comprise 
the scene. Next the computed size of the scene is printed follov^ed 
by several user supplied parameters which the program has utilized 
in carrying out its computations. The final listing is a summary 
of the v/ater body size distribution for the scenes (stratum) . 




LAKE AND POND STATISTICS 
SKYLAB/EREP S-192 DATA OF 12 JUNE 1973 
STRATUM: NORTH DAKOTA COTE AU 


S192 NCKTl- C4KCtA CJTA (<♦ . 2 £♦ 


Ci 5 GKT* 12*56+ 2.C35^ 12-06->72 


COCOCOCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCCCCCCCCCCCCCCOCGCCCOCCCCOaGCCCOCOaCCCCCCCCC 


TABULATICN 

CF REC0GM2EC WATER 

BCCIES 







AREA 

AREA 

LAT 

LCNG 

SCAN LINE 

. POINT 

(ACRES) 

IhECTARES) 

47.3C35 

99. 565 1 

285 

164 

1.227 

.497 

-47.2260 

59.6C6E 

287 

29 

6.136 

2.483 

47.3272 

59.5455 

291 

206 

47.861 

19,369 

47.21ie 

99.6C57 

292 

t 

1.227 

.497 

47.2403 

99.5936 

292 

56 

36.816 

14.099 

47.2113 

59,6C75 

294 

6 . 

1,227 ■■ 

.497 

47. 3286 

55.5C47 

297 

318 

1-227 

.497 

47 .2509 

59.5589 

298 

147 

44.179 

17.879 

47.2152 

99.5975 

302 

16 

19.635 

7.946 

47.2367 

99.584A 

3C2 

5^ 

1,227 

.,497 

47.2349 

99.5823 

306 

52 

9.818 

3.973 

47.2134 

99.5370 

307 

190 

1.227 

.497 

47.3377 

99.5222 

30 8 

233 

1.227 

,497 

47.3325 

99.5C88 

324 

230 

7.363 

2.980 

47.2012 

99.5799 

3 27 

1 

1.227 

.497 

47.2C07 

99.5781: 

329 

■ ■ 1 ■■ 

1.227 

.497 

47.2937 

99.5256 

329 

164 

29.453 

11.919 

47,2318 

99.5565 

2 23. 

57 

1.227 

,497 

47,2718 

99.5318 

335 

128 

142.355 

57.611 

47-2765 

99.5139 

345 

140 

1,227 

.497 

47.2760 

99.5171 

347 

140 

1.227 

.497 

47.3762 

99-4463 

361 

321 

2.454 

.993 

47,2097 

59.5382 

363 

30 

7.363 

2. 980 

47,3552 

99.4561 

263 ■ 

285 

1.227 

.49 7 

47.2422 

99 .5£9£ 

3 73 

91 

3.682 

1.490 

47.1957 

99. .5539 

375 

1C. 

7,363 

2. 980 

47.2436 

99,5068 

375 

94 

4.9C9 

1.937 

47-2065 

99.5196 

383 

32 

3.682 

1.490 

47.2872 

9 9.4 73 C 

384 

174 

4.909 

1.987 

47.2219 

99.5C68 

388 

61 

3.602 

1.490 

47.2,091 

99 .5120 

369 

: ■ 3 5 ■■■ 

2.454 

,593 

47.3614 

99.4259 

389 

306 

1.227 

.497 

47.3730 

99.4174 

391 

327 

30.680 

12.416 

47; 1855 

99.5110 

403 

3 

1.227 

.49 7 

47.2602 

99,4688 

403 

134 

1C3.C84 

41.718 

47.1915 

99.5065 

404 

. . 14' 

1.227 

.497 

47.1879 

99.5076 

4C5 

8 

6,136 

2.483 

47.1886 

S9.5C51 

407 

IC 

4. 909 

1.907 

47.3328 

99.42C6 

410 

264 

6,136 

2,4 6 3 

47.3665 

99,4016 

4 1 C 

323 

8.S9D 

3.477 

47.2591 

99.4366 

413 

. 206 ; . 

1.227 

.497 

47.2999 

99*4320 

417 

209 

1.227 

.49 7 

47,2857 

99.4349 

422 

186 

3,682 

1.49C 

47.2761 

99.4383 

424 

170 

6.136 ■ 

2.483 

47,3057 

99.4185 

427 

223 

2.454 

.993 

47.2086 

55.46 52 

435 

56 , 

12.272 

4.966 
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A7^2SIZ S?.A-!iCl A36 

A7.218S ?q,45AZ A AC 

47 ,281Z 5 9. 4L7C 442 

47.2517 99.4216 444 

47.3470 99.3722 448 

47.2415 99.4323 4-^9 

47.21CC 95..4.47C 452 ■ 

47.2516 99.5990 457 

47.1764 99.4557 462 

47.2612 99.4C6E 462 

47,2375 99.4179 465 

47.2847 99.3863 47C 

47.2794 ■' -2£c4 471 

47.2299 99.4143 473 

47.2597 99.3954 475 

47,2535 99 .39,68 477 

47.2821 ■ 9.9.3EC7 '■ v" ■ 477 

47.2042 99.42C6 461 

47.2234 99.4C36 427 ■ 

47.2551 99.2647 466 

47.1724 99.4211 498 

47.2452 99.273C 5C5 ' 

47.1773 95.4C92 527 

47,1932 99.4CC2 5C7 

47.2CC7 S9.396C 5C7 

47.3232 99.2237 SIC 

47.3289 99.22C4 51C 

47.2524 99.3621 511 

47.1557 99.2865 519 

47.3C71 99.3215 521 

47.1747 99.3933 524 

47.1943 99.2781 528 

47.1699 99.3785 53C 

47.1990 99.-724 53C 

47.1894 99.3767 532 

47.2539 99.34C3 532 

47.3018 99.3122 532 ■ /'. 

47.1753 99.2775 529 

47.1683 99.3EC5 .5*C 

47.2G06 99.3612 541 : 

47.2092 99.3513 646 

47.194 5 99,3575 54 E 

47,3046 99.29C2 552 

47.2142 99.34C2 554 

47.2100 99.2631 557 

47-ZG97 99,3276 559 

47.3284 99,27C6 559 

47.2146 59.333 6 - 56C 

47,3234 99.2714 561 

47.2142 99.3321 562 ' 

47,2041 99.Z6CZ 563 

47.2452 99.257C 563 

47.2297 99.3212 564 

47.285C 99.29CC 564 

47.2267 99.3152 56 6 .. 

47.3422 99.2556 566 

47.2496 99.3C69 567 

47.2232 99.2566 57C 

47.1762 99.3443 571 


131 


4.5C9 

1.987 

76 

'V 

2.454 

.993 

214 _ 


1.227 ■ 

.497 

106 


:■ 4.9C9 

1.987 

.135 


1.227 

.497 

3C5 


1.227 

,497 

119 


■■■ 6.136 - 

2.4 0 3 

- 65 


33.134 

L3.4C9 

21C 


309.253 

125.155 

ic 


28,223 

11.423 

159 


■ ■ 4.9C9 

1.98 7 

119 


1,227 

.497 

2C2 


■ 1.227 -j;. 

■■ .497 

194 


■ 0.590 . 

3.477 

ice 


25.771 

: . IC.43C 

161 


1.227 V 

.497 

151 


3.6B2 

' 1.490 

2 Cl 


■ -1.227 

.49 7 

66 


1 . 227 - ■■ 

■ .497 

1C2 


2.454 ■ 

.993 

158 


7.363 

: ' 2.9ac 

17 


1,227 

, .497 

147 


17.1B1 

■ .6^953 

29 


2CC.C33 

" 80 .953 

57 


17.181 

6.953 

7C 


31.907 

12.913 

206 


28.225 

' 11.423 

296 


2,454 

d-993 

164 


7.363 

2.98C 

0.6 


6.136 

2.4 03 

262 


56.451 

22.646 

31 


14.726 

.5.960 

67 


7.363 

2.980 

6C 


44.179 

17.879 

76 


16.408- 

7.45C 

6C ■ 


■ 11.C45 

■ • 4.47C 

173 


1.227 : 

■ .49 7 

257 


1.227 

.497 

38 ; 


1.227 

■ .49 7 

: 26 . 


■ 3.682 

1.49C 

83 


79.760 

32.282 

ICC 


39.270 

' 15.6,93 

. 75 


52.769 . 

21.356 

Z7C 


1.227 

.497 

U2 


31.907 

12.915 

281 


3.682 

1.49C 

1C6 


87.131 

35.262 

314 


3.682 

1.49C 

115 


2.454 

■ .993 

3C6 


'4.9C9 

1.987 

115 


2.454 

, .993 

273 


1.227 

.497 

34 5 


3.662 

1.490 

14 3 . 

■ ■; .. 

99.4C3 

40.228 

24C 


25.771 

LG, 4 30 

156 


,.126.401 ■■ 

51.155 

341 


2.454 

' .993 

179 


.. 11.C4 5 

4.470 

327 


1.227 

.497 


52 15.954 6-456 


. 
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■ ;:4'7,i6'04' ■ 

99.3467 . 

5.73 

39 

1.227 

.4 9 7 

47.2628 

■.■'99 .2010 

574 

24C 

6.136 

2.482 

■: ■4r.i£S5''-- 

99.3445 

575 

4C 

1.227 

.497 

. - 47.294C 

99.2737 

575 

26C . - , 

65.550 

26.305 

'-47.305ar' 

99.266C 

■■ '■" 576 

.. 281 

3.682 

■ 1.49C 

.47-3'l24\ 

99.2612 ■ 

577 

.293 ■ 

. . 11.945 : 

4.47C 

--V. 47-.242C;:- 

99.2CCC 

578 

■ 17C 

’ .227 

.497 

v' 47.3148 

99. 2978 

575 

■ 298 ■ ■ 

1,227 

.497 

47,3187 

99..252 5 

582 

■ 3C6 . 

1.227 

' .49 7 

.. 47.3204 , 

99,2.516 

: 502 ■ 

, 3C9 

: 2.4 54 

.593 

47.3347 

■ 99.2425 

502 

334 

2.454 

. . 5 53 

47.2722 

99.2747 

506 

. 226 

12.272 . ■ 

'■ 4.966 

. '47. 2796-: 

99.2593 

. 586 

274 ■/ ■ 

3,682 ■ . 

1.49C 

47.1567 

99.3389 

507 

24 

1.227 

. .497 

'47. '3259 ^ 

99.29C1 

" 508 

328 ■ ' 

4.905 

1.987 ■ ' 

'■ 47.2856 

99.264C 

585 ■ ■ ■ 

251 

13.495 

,5.463 

47.3182 

99.2456 

: 505 

3CE 

1,227 

.497 

47.1937 

99.315C 

59C 

■ 5E ' 

1.227 

.497 

'■ 47.2932 

99.2577 

591 

265 

1.227 

.497 

■ .1 47>1,870-; 

99.3167 

552 

75 ■' : 

107.993 

42.705 

47.33tfC 

99.2349 

■■■ ■ 593 ' 

23C 

1.227 

.497 

47.2978 

99.2521 

594 

274 

4.509 

1.987 

47.2804 

99.2563 

595 

258 

3.682 

1.49C 

47.2016 

9S.3C23 

598 

1C 7 

2.454 

.5 53 

47.3C70 

99.2418 

599 

252 

2.454 

.952 

47.2928 

99.2477 

6C1 

268 

34.361 

13.9C-6 

47.3042 

99.2413 

6C1 . 

ZEE 

6.136 - 

2.483 ■ 

■ : . : ■■ ■ 47.2730 

99.2563 

6C3 

9 234 

■ 1.227 

.457 

47.2815 

S9.252C 

6C3 

24 5 

■ 33.134 ■ 

13.4C5 ^ 

47.1826 

99.3C69 

6C4 

76 

19-635 

7.946 

47.2782 

99.2519 

6C5 

244 

11.045 ■ 

4,470 

47.2902 

99.2451 

6C5 

265 

2.454 

, .592 ■ , . 

; 47.1417 

99.3ZeC 

: 6C6 

5 ■ 

1.227 ■■ 

.497 

. 47.2061 

99.2885 

6C9 

119 

93.267 

37.74 5 

47.2298 

99.2741 

61C 

: 161 

4. 909 

1.987 

47.3X36 

99,2267 

6 1C 

3CE 

6.136 

2.482 

47.2601 

99.2549 

612 

215 . 

.. 88-3 58 ' 

35.758 

47.1635 

99.3C85 

613 

. 46 ,■ ■'■ 

4.9C9 

1.987 

47.2854 

99.2386 

614 

26C 

8.590 : 

3.477 

. 47,2481 

95.2587 

615 

195 . :: 

■■ 42.952 ■ 

17,383 

47.2995 

■99.2275 

617 

286 

■ : 2.454 . 

.592 

47.3207 

99.2156 

617 

323 

2.454 

• .592 

■ 47.3102 

99.2205 

. 615 

‘ 3C5' ■ 

'■ 1.227 ■■ 

.497 ■ 

47.2170 ’ 

99.2721 

615 

142 

9.B18 

3.973 

47.2392 

99.2596 

615 

181 

487,196 - 

197.168 

47-2757 

99.239C 

615 

24 5 

3.682 

1.49C 

47.1760 

99.2532 

621 

71 

2-454 

. .552 

47.2747 

95.2375: 

621 

:■■ 244 ■■ ■■■■ 

1.227 ^ .. 

.49 7 ■ 

47.2136 

95.2605 

624 

138 

1.227 : 

.49 7 

47.2456 

55.2488 

626 

195 

42.952 

17.383 

47.2673 

55.2365 

6Z6 

233 . - 

■ 1.227 

.457 

47-2315 

95.2547 

628 

171 

23.317 

9.436 

- 47.2628 

99.22 57 

6 28 

261 

7.363 : 

2 . sac 

■ 47.1919 ■ 

95.2761 

629 

■ ■ 1C2 

2.454 

.59 3 

47.2415 

9S.245C 

622 ■ 

' ISC 

: : 3.682 

1.49 C 

- 47.2660 

- 55.2311 

: 632 ■ ■ 

23 3 - 

1.227 

.497 

47.1374 

99.3C27 

6 23 

8 

1.227 

.497 

■ 47.2108 

99.Z6.C3 

'. 624 

137 

5 8.9C5 ■ ■■: / 

23.835 

■■ 47.1347 

99.3C22 

635 

4 

■ 6.136 

2.483 ^ 

■ 47.2237 - 

99.252C 

625 

■ - 16C 

6.136 

2.482 




^McyntfCBiLTiY of 
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<.7.2608 5S.231C 635 ' 225 . 3.682 I 


47.2813 

99.2194 

635 

2tl 

7.363 

2 

47.2816 

9 9 ,2182 

636 

V : 262 .. 

1.227 


47.2227 

99.2505 

637 

159 

. 6,136 . 

2 

47.2609 

99.2289 

637 

226 

■ ■■ 4.9C9 

: 1 

47.2196 

99.2512 

638 , , 

154 ■ 

■ ■ 2.454 , 


47,1376 

95.2955 

.640 . 

11 

1.227 


47.2826 

99.2115 

642 

266 

. . : 3.682 . 

. 1 

47.3113 

99.1933 

644 . ■ . 

317 

1.227 


47.2671 

99 .2059 

645 

■ 275 

1.227 


47.1457 

99.2837 

647 

26 ■ 

: 2.454 : 


47.2610 

99.2186 

647 

23C 

12.272 

4 

47.2545: . 

. 99.2213 

648 

219 

1,227 


47.2778 

- 99.2CS1 ■ 

648 

260 • ■■ 

2.454 


47.2634 

99.2152 

649 

235 . 

1.227 


47.2736 

99.2094 ; 

649 

253 

33.134 

13 

47.2609 

99.2i56 

65C 

231 

1.227 


47,2660 

99.2127 

650 

■■■■■ 24C 

■ ■■ ■ 47.861 

19 

47.1557 

99.274C 

^ ■ 651 ■ 


4,909 ■ ■ : 

1 

47.2549 

99.2179 

651 

221 

2.454 . 


47.1974 

99.2454 

656 

: 122 

1.227 


47.1983 

99,24 38 

657 

124 

1.227 


47.2637 

99.2039 

658 

■ 239 ■ 

1.227 


47.2694 

99.2026 

658 

:■ 245 

1.227 


47.2950 

99.1881 

658 

294 

1,227 


47,2995 

99,1825 

661 

303 

■V . 2.454 


47.1515 

99.2651 

662 

: 44 

2.454 


47-1681 

99.2557 

. 662 : 

■ - • 73 

2.454 


47.2605 

99,2035 

■ 662 

235 

6.136 

2 

47.2719 

99.1971 

662 

255 

2.4 54 


47.2534 

99.2065 

663 

223 : 

. 29.453 ■ ■ 

.11 

47.1534 

95.2620 

664 

•48 

£.590 

3 

47.2809 

99.1889 

665 • ■ 

.. ■ 272 

■ 1,227 - 


47.2571 

99.2CC3 

66 7 

231 

. 9.aia :■ 

- 3 

47.2C90 

99.2265 

6 68 

147 

76.006 

• 30 

47.2475 

^ 99.2037 

669 

215 

9.618 

. . 3 

47.1505 

99.2555 

672 

46 

4.9C9 

, 1 

47.1814 

95.2359 

674 

ICl . . , 

2C.862 : 

... . 8 

47.1946 

99.2264 

676 

' 125. 

1,227 ■■ -V 


47.2677 

99-1652 

676 

253 

11.045; ■ 

4 

47.1260 

99,2642 

677 

C 

4.9C9 

1 

47.2806 

59.1769 

677 

276 

1.227 


47.1719 

99.2372 

678 

86 . 

. 67.496 

27 

47.2402 

■ 99.1976 ■ 

■■■ 679 ■ 

206 ■ 

6,590 


47,2034 

99.214 3 

683 

143 

11.045 

4 

47.2576 

99.1837 

683 

226 - : 

7.363 

if 

47.2366 

99.1944 

684 

2C2- 

4.509 

1 

47.1305 

99.2535 

685 ; 

16 

■- 3.682 

1 

47.2C35 

99.2122 

685 

144 

1.227 


47.2143 

99.2061 

■ 6B5 

163 

95.721 

. .36 

47.1947 

95.2161 

686 

, 129 

3.682 

1 

47.2415 

99,1897 

686 

211 

16.408 

7 

47.2494 

59,1822 

689 

226 

1-227 


47.2549 

59,1781 

65C 

236 

2.454 


47.2076 

99.2027 

692 

154 

; 65.041 

■26 

47.1685 

99.2207 

696 

. 87 : ■ r: 

^ ■ 1.227 


47.1948 

99.2C59 

696 

123 , 

. 1.227 


47.2260 

99.1856 

697 

193 

, 2.454 


47.2402 

59.1747 

697 

227 ^ 

1.227 



<^7.1359 99.236C 699 31 14.726 


.99G 

,98C 

.497 

.483 

.987 

.993 

.4 97 

.49C 

.497 

.497 

.953 

.966 

.497 

.993 

.497 

.409 

.49 7 

.369 

.987 

.9 9 3;: 

.497 

.497 

.497 

.497 

.497 

.993 

.993 

.9 93 

.483 

.99 3 

.919 

.477 

.497 

.97 3 

.792 

,973 

.987 

.44 3 

.497 

.470 

.987 

.4 9 7 

,316 

,4 7 7 

.478 

.960 

,987 

.490 

.497 

.738 

-49C 

.4SC 

.4 97 

.99 3 

,3 22 

.497 

.497 

,993 

.497 

.96C 



^7.2049 

59.1971 

699 •' 


152 


7.363 

2.980 

^7.2<i96 

99,1647 

706 


' ' 233 ' '■ 


1.227 

. ' .497 

A7.Z3^6 

59.1721 

707 


207 


1.227 

,497 

A7.1379 

99.2257 

- 7 OS . ' ■ 


26 - 


28,225 

; 11.422 

A7.1302 

99.229C : 

709 


25 


44.179 

17.875 ':■■■• 

AT. 2207 

59.1723 

■ 7C9 


201 


4.5C9 

1.987 

A7.1A6& 

99.2187 

71C 


54 . 


1.227 ■ : 

.497 

A7.178t 

99.2C07 

710 


lie 


9.818 

3.572 

■,: A7^^2182 :;: 

99.169U 

' 719 .'r 


182 


- 3.682 ..- 

■ 1.490 

■ A7.I651 

.99.1960 ■ 

7 22 


91 ■ 


1.227 

.49 7 - 

- A7.2i2i::: 

99.1675 

724 


174 


: ■ 6.590 ^ : 

■E.477 : 

A7.2A12 

99. 14 9 C , 

^ 726 


226 


2,454 . 

- .992 

47.2228 

99. 1C 29 

. 726 


369 


2.4 54 - ■ .'■ 

.992 

47.3250 

99.1C16 

.726 


373 


6.136 - 

, : 2 .4 8-2 

47.1172 

99.Z18C 

727 


- 9- • 


2.454 

'■ .'5 9 2 ' . 

A7.1395 - 

•5 9. '15 51 

73 ? . ‘ 


■ ‘.52 


34.361 

■ 13 .■90 6 . 

47.1771 

; 9 9.1739 

727 ■ 


HE 


-■ .'■ 11.04 5 

4-470 ' ' 

47.1159 

99,2G74 

728 , 

T' - 

. 11 


4.9C9 

1*9.8/ 

47.2526 

99.1292 

739 


251 


11.045 

' : 4- . 7 ?c . 

V 47.1768 

■99 .1699 ' 

■7'40'. 


122 


1,227 

>7 

47.1791 

99.1687 

.741 


122 


1.227 

.4-7 

A7.188i 

99.1622 

741 : 


14C 


1.227 

,49 r . 

47.1655 

99.1742 : 

742 


ICC 


2.454 

.9 9 2 

47,1767 

99,1669 

: 742 ' 


- 122 


'■' 2.454 : ■ 

.592 ■■': 

47.1912 

99.1577 

745 


146 


: 1,227 

.497 '■ 

A7.14Ce 

59.1821 

749 


59 


6.136 

2.463 

L 47.2C52. 

: 99.1457 

■■■ 749 


■ 172 


2.454 

' '- .992 

47.1451 

95.1756 

753 


: 68 


6.136 

2. 482 ■ 

47.1112 

99*1537 

754 

.... 

. 9 


31.9C7 . 

12.912 

47.146G 

99.1741 

754 V 


: 70 ■ 


2.454 

.992 ■ 

47.2098 

99.125C 

7 57 


183 


2.454 

.592 , 

,47.1256 

59.1784 

761 


■ 37 ■ 


2,454 : 

-592 

|47.1212 

59.1785.. 

; 762 


30 


69.550 . ' 

28-305 

47.2145 

59 .1251 

764 


194 


4 . 509 

1.987 

,4 7. 142 7 
47.1365 

99.1606 

769 


70 


1.227 ' 

,457 ■■ 

99.1529 

779 


63 


: 23.317 

• 9.436 

47.1567 . 

99*1282: ' 

■ 7 82 


ICC 


12.272 

4.966 ■■ 

47.15CG 

5 9. 1401 ' 

785 


■ 89 ; 


6.136 

2.483 

: 47.1234 

59.1510' 

7 69 


44 


115.G3B 

46.175 

47.2689 

55*0627 

794 


: ■ 3C1 


£.590 

"V 2.477 

47.1472 

95.1214 : 

795 


88 


1.227 

.497 

- 47. 1460 

99.1251 

■^:'79a 


67 


1.227 . 

.497 : 

47.1078 

99.1455 ' 

. 402 . . 


: 22 


3G.6B0 

12-416' 

47.1602 

59.1159 

802 : 


114 


‘ ^ 15.954 

6.456 

"■47.126 8 

.59.1261. 

'.806. 


74 


44.179 

17.875 

47.1737 

95.1021 

8C9 


14C 


. 1.227 

; .4 9 7 

. 47.1747 ; 

; 5 9-.. ICC 6- 

810 . 


142 


4.909 

1,98 7 

47.1257 

99.1261 

812 - 


57 


34.361 

13.906 ■ 

■■47.1C37 

99.1215 

819 


21 


1.227 

.497 . 

47-1468 

99*1.061 ; 

820 


97 


1.227 

.497 

.47.0929 

99. 13-55 

*. 821 


3 


■ 18.408 

7.450 ■ 

”47 . 1290 

: 99.1079.. 

' -.'828 


69 


2.454 

.593 

AT. 1655 

99.0873 

.828 


133 


1,227 

.49 7 '■ ■ 

47.1237 

.99.1059 

■ 829 


6C 


2.454 

.553 

47.1247 

99.1062 

•: 822 


■ 63 ' 


1-227 

■ .497 

47.1C17 

59.1182 

82:2 


23 


85.904 

34.765 

47.1045 

99 . 1 16 6 

822 


. 28 ... 


- 2ii4 54 

- .993 

47,1226 

■' 99. 1054 

824 " 


- 6C 


■ ■■ 2.454 ■ ■ 

■-■■r .553 

. 47.1643 

99.0768 

S29 


135 


1.22 7 

.45 7 

47-1669 

99,C74T 

840 


14 C 


1.227 

-49 7 




PGOE 




A7.1670 

.9^.0721 

84 2 


141 

1.227 

.457 

«7.1312 

•59.C9C3 

844 


79 

y. '- > ■■/ 1.227 

.497 

47.1036 

9 9. 1C 4 9 

■'84 5 


31 

■■ - yy 2,454' ■ 

,993 

A7.117S 

9 9. C 96 6 

845 . 


56 

.15,954 

6,4.56 

^7.1032 

95>IC31 

. 647. ■ 


'■■V 

.■.4,sc9 - 

1.967 

A7.ll6a 

99.0953 

84 7 


C'C 

8.590 ■ 

3.477 

47.1366 

99.0622 

-84 8 


90 

; ye.59G -y- 

3,477 

47.U1E 

99.C961 

849 


4 7 

3.682 

1.490 

4TU43S 

99.CT6G 

851 


■ ilC4 

2.454 

,953 

47.1627 

9 9.0653 

851 


137 

3.682 

■y/ ;- :i-.4 9C; "■ 

47.1021 

99.CS66 

852 - ■ 


31 

2.454 

.993 

47.1395 

^9.C764 

853 


57 

. , 11,04 5 

. 4.4 70 

47.1501 

99,0654 

■ 8 54 


116 

7.363 

2.5&C 

47,0910 

99,0967 

6 58 


14 

6.590 

3.477,: 

47.1047 

9 9.C9.C9 

- 35 6.'-:' 


^'38 ■ 'y 

■ y ■ y y ■.6,590 " 

■ 3,477 ■ 

47.1173 

99.0639 

. 858 


. ■ 6C 

2.454 

.993 

47,1412 

: 9 9.3703 

. 856 ■ 


■ ICZ 

3.682 

1-4 5 C . 

47,1231 

99.C765: 

66C 


71 

;■ 4.9C9 

1.93 7 ■' 

47.1446 

99.C6S4 

■ 661 


IC9 

1.227 

. .497 

47,1180 

9 9... 0763 . 

£6 5 


■■ ■■.y- .64... 

-.1,227 

■ y : ■ *497 ■ 

47.1385 

99.G647 

865 


ICC 

1,227 

.497 

47il437 

9 9, C 6 16 

6 65 


IC9 

7,363 

2*98 C 

4 7.1 GO 2 

99.C632 

■ ■66 8 ■ 


34 

3.682 

•.■l,490"y 

47. os's 2 

99.G907 ' 

■ 669 


8 

3.682 

• 1.490 

47.1240 

.99,0686 

669 


76 ■■ 

5.818 

3.973 

47.0976 

59.CE15 

871 


21 

2.454 

.993'' 

47.CS85 

9 9,C65e 

87 2 


^ 15‘ 

14.726 

y- 5.960 

47.136E 

95.0575 

872 : 


ICC 

17-181 

6.953 

47.0959 

99. 0765 

877 


3C ■ 

E.550 

3.47 7 

47.1153 ■ 

99,0655 

877 


'■ 

-y^- - y y: : 1.227. 

y^ ,497 y 

47.1166." 

99.C626 

677 


; 7C 

4.909 

1.987 

47.1197 

5S.062C 

67 6 


■ ■■ .'72 y 

2.454 

: .993 . 

47.0865 

99.G795 

B79 ; 


: ' y 15 . . 

23.317 

9.436 

47.1061 

9 9, C 6 7,6 

660 


49 

1.227 

.497 

47.1116, 

99,0644 

y 6EC-, 


yyy ;^ 59 v- y 

. • y 2. 4 54-, 

y.y ■ :y .99 3- 

47.0945 : 

9 5,0722 

6 El 


" ■29'.'-- 

e.,59C‘ 

.,■■ ■ 3.477 

47.0979 

99.G712 

661 


■ 3 5 

45 .406. 

18.376 

47.2331 

9 E. 9949 

6 El 


. 272 

14.726 

5,960 

47.1CQ5 

9 5, 06 a 7 

6 62 


4C 

1.227 

■ .497 ■ 

47i‘L380 

■95:,C46t 

y 66 3 ' 


:. . 1C6 . 

y y y y.."6,:454 

.99 3 

47.1283 

99.0469 

; : : £68 ■ ■ 


91 

7,363 

2.98 0 

47,1347 

99.04 12 ■ 

/ 890 .. 


1C3 

6., 136 

2.483 

47.0966 

99.3616 

651 


y ■ 37 

-3-682;. 

■ 1.490 

47.0989 

99.0554 

892 


41 

3.682 

, 1.4 90 

47.1248.^' 

9 9,0407 

896 ■ ^ 


■■yy £E„- 

1.227 

.497 

47.Q76C 

■9 9.0 642 

9CC 


„ . A- 

■ : 29.453 

■ 11.9:19 

47.1186 

99.C43C 

90C 


75 

11,045 

y : 4.47C 

47.0944 

59.0517 

902 


■ 37 y 

2.454 ■ " 


47.0740 

99.C612 

9C4 


2 

3.682 

1.49C 

47.083E 

99,0536 

9C6 


;,y, -2C.... y.- 

2,454 

■ .9^z: 

.47.0831 

99.0376 

■ 922 


..^25. - 

y- yyyy ■' 4G,497 . 

■ ■ i6,3 8'9y 

47.0715 

99.0432 

522 


. .5 •• 

19.635 

7,946 


SCENE N.C* COTEAU LINES 285 THRU A25 PCLMS I TFRU 28C 

y LINES ^26 THRL 530 PCINTS 1 THRU 2A1 

LINES 531 THRU 635 PCINTS 1 ThBU 3ZA 



10:2 


L INES 

636 

THRL 

740 PCINTS 

1 

TPBU 

272 

LINES 

741 

ThRL 

845 PCINTS 

I 

TPBU 

2C1- ■ 

LINES 

8A6 

THPL 

S24 PCINTS 

I 

TPPU 


scene: 

AREAf= 

. 34 8 

: SC. ''r &■■■■■ 





_ * • 

<5C2 

sc. KP. 




PIXEL 

LENGTH 

!= :-Ti 

.S44 PET EPS PIXEL 

HlCTh- 

69.C2L PETERS 


KCGE= ■ ^CR^'AL ■ ' ■ 

PCIMS CCLNTEC IF VCLTAGE IS GREATER Tl-AN CP ECUAL TC C AAC 
LESS THAN CR ECLAL 7C 1C ’ 
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S192 NCRTf- CAKCTfl DATA »<* 2 S* C* 9 Gf^T* 12=^56^ 3.025A 12-06-72 

OCOOOOOOCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCGeCCCOCCCCCCCCCCCCCCCC 


SCENE 

N.C. CGTEAf 

Lines 

26 b 

THRU 

925 

PUNTS 

1 

THRU 

sec 



LINES 

A26 

IHRL 

53C 

PUNTS 

1 

THRU 

391 



LINES 

531 

THRL 

6J5 

PUNTS 

1 

THRU 

229 



LINES 

636 

ThRL 

790 

PUNTS 

1 

THRU 

273 



LINES 

7Al 

IHRL 

E95 

PUNTS 

1 

THRU 

201 



LINES 

S46 

ThRL 

529 

PUNTS 

1 

THRC 

1C6 


CISTRIELTION CF RECCGNIzeU t>ATER BCCIES IN THE SCENE 
PY AREA 


AREA 

(ACRES) 

AREA 

(HECTARES) 

FRECLEKCY 

.25 

TC 

.5C 

.1C 

TC 

.2C 

C 

.50 

TC 

l.CQ 

.20 

TL 

.9C 

0 

1-00 

TC 

2. CO 

.9C 

TC 

.61 

58 

2.00 

TC 

3.CC 

.61 

TC 

1.21 

51 

3. CO 

TQ 

9.CC 

1.21 

rc 

1.62 

27 

9.00 

TC 

6. CO 

1.62 

TC 

2.93 

23 

6.00 

TC 

6.CC 

2.93 

TC 

3.29 

31 

B.CQ 

TC 

it.cc 

2.29 

TC 

9.C5 

19 

10. CO 

TD 

15. CC 

i.C5 

TC 

6.C7 

19 

15.00 

TC 

2C.C0 

6.07 

TC 

E.C9 

12 

20. CC 

TC 

25.CC 

6.C5 

TC 

1C. 12 

9 

25.00 

TC 

2C.CC 

1C. 12 

TC 

12.19 

a 

30.00 

TC 

9C.CC 

12.19 

IC 

16.19 

13 

90.00 

TC 

5C.CC 

16.19 

TC 

2C.23 

IG 

50.00 

TC 

75. ce 

2C.23 

TC 

3C.35 

7 

75.00 

TC 

ICC.CG 

2C.35 

TC 

9Cc97 

8 

icQ.ac 

TO 

15C.CC 

9C.97 

TC 

6C.7C 

5 

150.00 

TC 

2CC.CC 

6C.7C 

TC 

EC. 99 

C 

CVER 2C0.CC 

CVER 

EC.' 

99 

3 


i. 

t 
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LAKE AMD POND STATISTICS 
SKYLAB/EREP S-192 DATA OF 12 JUNE 1973 
STRATUM; NORTH DAKOTA DRIFT PLAIN (WEST) 


SI92 NCRTh CAKCTA DiTA 2 0# 9 12*56* 3.G354 12-06-72 

CCCCCCCOCGCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCCCCCCGCCCCCCCCCC 


TAEULATICA CF RECOGNIZED N4THR RCCIE5 


AREA AREA 


LAT 

LCNG 

SCAN LINE 

PCINT 

(ACRES) 

(HECTARES) 

^7.AS35 

S9.395C 

348 

5C4 

4,909 

1.987 

47.A657 

99.4C59 

351 

474 

. 2.454 

.953 

47.4761 

99.398C 

353 

493 

29,453 

11.519 

47.4702 

99.4CC3 

354 

483 

1 .227 

.497 

47.4688 

95 ,4CC1 

355 

481 

1.227 

.497 

47.5319 

99.2E97 

4 28 

62C 

2.454 

.993 

47.5347 

99.2800 

436 

628 

6.136 

2.483 

47.4592 

59.2582 

499 

52C 

1.227 

.497 

47.4483 

95.251C 

512 

5C6 

3C.680 

12.416 

47.5C32 

99.1627 

568 

624 

1.227 

.497 

47.5028 

99.16C9 

S7C 

624 

1.227 

.497 

47,3817 

95.1627 

635 ; 

437 

22.C90 

a,94C 

47.4623 

99.I1CC 

642 

581 

19,635 

7.946 

47.3702 

99.1565 

643 

434 

51.542 

2C.659 

47.3594 

99.1254 

674 

413 

13.499 

5.463 

47.3814 

99.1168 

6 8C 

4 54 

117.811 

47,678 

47.2335 

95 .0886 

734 

391 

15,564 

6.456 

47 .4396 

95.C1C2 

752 

564 

1.227 

.497 

47.4114 

98 .9668 

eic 

557 

1,227 

.497 

47.3119 

9E.9554 

837 

393 

1.227 

.497 

47.3913 

9E.9434 

844 

535 

36C.795 

146. C14 

47.4484 

98.8507 

9C3 

668 

6.136 

2 .483 


scene 

-CRIFT PLAIN WEST 

LINFS 

285 THRL 

425 

PCUTS 

381 

THRU 

67 



LINES 

426 THRL 

530 

PCUTS 

342 

THRU 

67 



LINES 

531 THRL 

635 

PCUTS 

325 

IHRU 

67 



LUES 

636 THRL 

74C 

PCUTS 

274 

THRU 

67 



LINES 

741 THRL 

34 5 

PCUTS 

2C2 

THRU 

67 



LUES 

846 THRL 

924 

PCUTS 

107 

THRU 

67 


SCENE AREA= 

480 

SC. RI. 







= 

1244 

SC. Kf*. 







PIXEL L£NGTHi= 71 

.944 PETERS 

PIXEL 

HICTH= 

69.C31 PETERS 


RCDE^ NCRNAi 

PCIMS CCLMLt. IF VCUAGE IS GREATER Tt-AN Cfi ECUAL TC C A^D 
LESS THAN CR ECLAL TC 1C 


in in %n in kn \n 
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S192 KCBTh CAkCT/S C^Ta W* 2 S* c* 9 GM* 3.C25<> 12-36-72 

ccoccacccccccccccccccccGcccccccccccccccccccccccccccccccacocccoccccccocccccccccc 


SCENE 

CBIFT PLAIN WEST LINES 

265 

THRL 

425 

PCINTS 

581 

ThRU 

675 


LINES 

A26 

Thrl 

53u 

PC IMS 

• 342 

ThfiU 

675 


LINES 

531 

ThKL 

6i5 

FCIMS . 

325 

TERU 

675 


LINES 

6 36 

TNRL 

74C 

PL IMS 

274 

TWPU 

675 


LINES 

741 

ThKL 

d45 

PCINTS 

2C2 

TEBU 

. 675 


LINES 

346 

Ih-iL 

924 

FCIMS 

107 

TEHU 

675 


CISTRIBLTION CF BECCGMZEC V*aTER liCClES IN Thb SCENE 
BY AREA 


AREA 

(ACRES } 

AREA 

(HECTARES) 

FRECLENCY 

.25 

TC 

.5C 

.1C 

TC 

.2C 

0 

.50 

TC 

l.CQ 

.2C 

TC 

• 4C 

0 

l.OC 

TG 

2.CC 

.4C 

TG 

.El 

8 

2.00 

TO 

3.CC 

.El 

TC 

1.21 

2 

3. CO 

TC 

4.CC 

1.21 

TC 

1.62 

C 

4. CO 

TC 

6.CC 

1 .62 

TC 

2.43 

1 

6.CC 

TC 

6.CC 

2.43 

TC 

2.24 

2 

8. CO 

TC 

1 L * U C 

2.24 

IC 

4. (.5 

0 

10. CO 

TC 

15. CC 

4.C5 

TC 

£ .C7 

1 

15.00 

TC 

2C.CC 

6 .C7 

TC 

£.'9 

2 

2C.0C 

TO 

25. CC 

£.C5 

TC 

1C.12 

L 

25.00 

TC 

3C.CC 

IC.12 

TC 

12.14 

1 

30. CO 

TC 

4C .CC 

.12.14 

TC 

16.19 

1 

40.00 

TC 

5C .CC 

U.19 

TC 

2C.23 

0 

5C.C0 

TC 

75. CC 

2C.23 

TC 

3C.35 

1 

75.00 

TC 

ICC .CC 

20.35 

TC 

4C.47 

C 

100. CO 

TC 

15C.CC 

4C.47 

TC 

6C.7C 

1 

X5c,ca 

TC 

2CC.CC 

ec,7c 

TC 

EC. 94 

Q 

CVER 2C0.CC 

CVER 1 

EC. 94 

1 


\n vn m in in in 


LAKE AND POND STATISTICS 
SKYLAB/EREP S-192 DATA OF 12 JUNE 1973 
STRATUM: NORTH DAKOTA DRIFT PLAIN (EAST) 


S19Z MORTl- CAKCTA CiT/i 2 S* C# 9 Gf*!* 12*56* 3.0354 12-06-73 

CCGOODCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCOOCCCCCOOCCCCCCCOC 

TjSEULATICN CF RECOGMZED hATER BCCIES 


AREA AREA 


LAT 

LCNG 

scan lime 

FCIM 

(ACRES) 

(HECTARES) 

47-0705 

59.C321 

93C 

20 

1.227 

.497 

47.0691 

99.0333 

934 

e 

1.227 

.497 

47.4120 

98.8248 

946 

613 

49.C88 

19.866 

47.3995 

98.8213 

949 

59C 

3.682 

1.490 

47.4195 

96.8149 

954 

627 

3.682 

1.490 

47.2795 

96.G53C 

955 

362 

11.045 

4.470 

47.1659 

96.9541 

950 

164 

1.227 

.497 

47.1671 

96.9513 

96C 

187 

1.227 

.497 

47.4427 

96.7957 

96C 

670 

3.682 

1 .490 

47.3253 

96.8548 

967 

467 

4.9C9 

1.987 

47.0743 

98.9832 

96C 

32 

1.227 

.497 

47.3297 

96.6329 

966 

4 62 

6.136 

2.483 

47.2600 

98.8692 

969 

361 

1.227 

.497 

47.2507 

96.3668 

992 

36C 

1.22 V 

.497 

47.2616 

96.0652 

992 

365 

1.227 

.497 

47.2577 

96.6653 

994 

359 

3.682 

1 .490 

47.2634 

98.8621 

994 

369 

1.227 

.497 

47.4170 

96.7688 

. ICCC 

642 

3,682 

1,490 

47.3235 

96.8143 

1CC2 

495 

2.454 

.99 2 

47.3347 

98.8137 

1CC2 

497 

1.227 

,497 

47.1706 

98.9C22 

1C06 

211 

12.272 

4 .966 

47.3121 

96.0172 

ICIl 

461 

15.954 

• 6.456 

47.2736 

96.6349 

IC15 

395 

1567.128 

634.217 

47.3428 

. 96.7936 

1017 

517 

4.9C9 

1.987 

47.2683 

96.8317 

1021 

366 

3.682 

1.49C 

47.2994 

96.8131 

1022 

443 

1.227 

.497 

47.1107 

96.9176 

1C24 

113 

1.227 

.497 

47.2762 

96.8242 

1CZ4 

403 

1.227 

.497 

47.2773 

96.6236 

1024 

4C5 

1.227 

.497 

47.1446 

96.8914 

1031 

175 

3.682 

1 .49C 

47.2278 

96 .8442 

1031 

321 

2.454 

,993 

47*2259 

96.0444 

1032 

318 

181.625 

73.504 

47.2628 

90,8226 

1033 

363 

6.136 

2.483 

47.2C09 

98.8565 

1034 

275 

1.227 

.497 

47.2409 

96.8339 

1034 

345 

26.225 

11.423 

47.2002 

96.8117 

1034 

414 

1,227 

,497 

47.1988 

96.6556 

1C36 

272 

6.136 

2,483 

47,2242 

96,8257 

1036 

334 

31.907 

12.912 

47.2427 

96.63C8 

IC36 

349 

3.682 

1.49C 

47,1704 

96,6696 

1038 

223 

4.9C9 

1,987 

47.2195 

96.8316 

1048 

313 

1.227 

.497 

47.2421 

96.6179 

1049 

353 

23,317 

9,436 

47.1637 

98.8611 

1C5C 

216 

1,227 

.497 

47.2251 

96 ,6254 

IC51 

324 

4. 909 

1 .987 

47.2310 

96,eC16 

1071 

342 

1.227 

,497 

47,2218 

9E.76C2 

1097 

336 

1,227 

.497 
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AT -220 3 

98.7769 

lici 

315 

6,136 

2.483 

47.1604 

98.7678 

Ills 

251 

765.769 

309.907 

47.0046 

96.6U8 

1137 

111 

l.??7 

.497 

47.1363 

98.7313 

1192 

22 3 

31.907 

12.913 

47.2292 

98.6757 

1195 
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